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Abstract 

The subject matter of this paper concerns the derivation of the finite Larmor 
radius approximation, when colhsions are taken into account. Several studies 
are performed, corresponding to different colhsion kernels. The main motivation 
consists in computing the gyroaverage of the Fokker-Planck-Landau operator, 
which plays a major role in plasma physics. We show that the new collision 
operator enjoys the usual physical properties ; the averaged kernel balances the 
mass, momentum, kinetic energy and dissipates the entropy. 
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1 Introduction 

Many studies in plasma physics concern the energy production through thermonu- 
clear fusion. In particular this reaction can be achieved by magnetic confinement i.e., 
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a tokamak plasma is controlled by applying a strong magnetic field. Large magnetic 
fields induce high cyclotronic frequencies corresponding to the fast particle dynamics 
around the magnetic lines. We concentrate on the linear problem, by neglecting the 
self-consistent electro-magnetic field. The external electro-magnetic field is supposed 
to be a given smooth field 

E = -V,<P, B' = ^h{x), \h\ = l 

when e > is a small parameter, destinated to converge to 0, in order to describe strong 
magnetic fields. The scalar function cf) stands for the electric potential, B{x) > is 
the rescaled magnitude of the magnetic field and b{x) denotes its direction. As usual, 
we appeal to the kinetic description for studying the evolution of the plasma. The 
notation = /^(t, x,v) > stands for the presence density of a population of charged 
particles with mass m and charge q. This density satisfies 

dtf + V ■ V.r + -iE + vAB')- V.r = Qif), (t, X, v)eR+xW'x (1) 
m 

/"(O, X, v) = f f ), (x, v)eR^ X M3 (2) 

where Q denotes a collision kernel. The interpretation of the density is straight- 
forward : the number of charged particles contained at time t inside the infinitesimal 
volume dxdv around the point (x, v) of the position-velocity phase space is given by 
/^(t, X, f )dxdf . The equation ([T]) accounts for the fluctuation of the density due to 
the transport but also to the collisions. We analyze here the linear relaxation operator, 
but also the bilinear Fokker-Planck-Landau operator. 

When neglecting the collisions the limit model as e \ comes by averaging with 
respect to the fast cyclotronic motion [IS], [211 El [D El El IH IS] • The problem reduces to 
homogenization analysis and can be solved using the notion of two-scale convergence 

[MlEllls]. 

We point out that a linearized and gyroaveraged collision operator has been written 
in [22], but the implementation of this operator seems very hard. We refer to [E'j'PJ for a 
general guiding-center bilinear Fokker-Planck collision operator. Another difficulty lies 
in the relaxation of the distribution function towards a local Maxwellian equilibrium. 
Most of the available model operators, in particular those which are linearized near a 
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Maxwellian, are missing this property. Very recently a set of model collision operators 
has been obtained in [TS], based on entropy variational principles [7]. 

We study here the finite Larmor radius scaling i.e., the typical perpendicular spatial 
length is of the same order as the Larmor radius and the parallel spatial length is much 
larger. We assume that the magnetic field is homogeneous and stationary 

B'= fo,0,- 



e 

for some constant > and therefore ([T]) becomes 

where Uc = qB/m stands for the rescaled cyclotronic frequency. The density is 
decomposed into a dominant density / and fiuctuations of orders e^e^, ... 

r = f + ef + e'f + ... (4) 

Combining ([3]), (H]) yields, with the notations x = {xi,X2),v = {vi,V2), '^v = {v2, ~Vi) 

Tf:=v-V^f + oo,^v-V^f = (5) 

dtf + v,d,J + ■ + rf = Qif) (6) 
m 

The equation ([5]) appears as a divergence constraint 

div,,,{/(U,0,u;,^T7,0)} = 0. 

Equivalently, ([5]) says that at any time t the density /(t, - , ■) remains constant along 
the fiow associated to v ■ V^- + Uc'^v ■ 

dx - dXs ^ dF ^- m . 

and therefore, at any time t, the density /(t, ■, ■) depends only on the invariants of ([7]) 

f{t,x,v) = g [t,xi^ ,X2 ,X3,r = \v\,V3 

The time evolution for / comes by eliminating in ([6]). For doing that, we project onto 
the kernel of T, which is orthogonal to the range of T. In order to get a explicit model 



for / we need a simpler representation for the orthogonal projection on ker T. Actually 
this projection appears as the average along the characteristic flow ([7j). Denoting by 
(■) this projection, we obtain 

(dtf + v;d^J + ■ Vj) = {Qif)) , {t, X, t;) G M+ X R3 X Ml (g) 
\ ml 

By one hand, averaging dt + v^dx^ + -^E - Vv leads to another transport operator 

dtf + v^dxj + ^E ■ Vj) = dtf + ■ V^/ + vsdxj + ^ {Es) d,,f. 

m / B m 

The key point here is to choose as new coordinates the invariants of ([7]) and to ob- 
serve that the partial derivatives with respect to these invariants commute with the 
average operator. More generally, for any smooth vector field ^ = {C,x,C,v), we obtain 
the following commutation formula between the divergence and average operators, cf. 
Proposition 13.31 

(div.,,0 = div^{(^^+^) + (^^- ^) ^ - (^^- ^) ^] +dx,{U 

I \ CUc / \ \V\ / (^r\v\ \ If I / Ur\v\ ' 

' 1^1 / 1^1 

By the other hand we need to compute the average of the collision kernel Q which is 
a more complicated task. It is convenient to focus first on the relaxation Boltzmann 
operator [21] 





Uc\v\ 


V 




\v\ 





QB{f{t,x,-)){v) = - I s{v,v'){M{v)f{t,x,v') - M{v')f{t,x,v)} dv' 



T 



where r > is the relaxation time, s{v, v') is the scattering cross section and M is the 
Maxwellian equilibrium with temperature 6 

We need to average functions like {x,v) — t- J^3C{v,v')f{x,v') dv', where C{v,v') is 
a given function. Since the invariants of the flow {X, V) combines x and v, we get a 
position- velocity integral operator cf. Proposition 14.21 

(/ C{v,v')f{x,v')dv'){x,v) = ul / C{\v\,V3,\v'\,v'2, z)f {x',X3,v') dv'dx[dx2 
\Jr3 I Jr^J^s 
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with z = UcX + -^v — {(jJcX' + -^v'). We prove that averaging Qb will lead to a position- 
velocity integral operator of the same form 



{QB)f{x,v): = {QB{f)){x,v) 

= / 5(H,t;3,PU^,^){M(t;)/(^,X3,t;')-MK)/(x,t;)}dt;'dx;dx^ 



r 



(see Theorem 11.11 for the definition of S). Observe that {Qb) is global in but 
remains local in X3. In particular it satisfies only a global mass balance, which comes 



easily by Fubini theorem and the symmetry of S cf. Remark 14.61 

{Qb) f{x, v) dvdx = 0. 




In the case of the relaxation operator Qb obtain the limit model 

Theorem 1.1 Assume that the scattering cross section satisfies (1221) . fl271) and that 
E{x) = -V^(j){x),(f) e W^^^im:-^). Let us consider > 0, f G L\R^ x M^) n 
(M~^dxdt') and denote by the weak solution 0/ (j3]), (j2]) with Q = Qb for any 
e > 0. We assume that {f'^)s>o is bounded in (M^-'^dxdf )) . Then the family 

(/^)e>o converges weakly in L°° (M+,L^ (M~^da;df)) to the weak solution of 

dtf + ^^-V^f + vsd^J + ^{E,)d,J={QB)f, {t,x,v)eR+xR^xW' (9) 
B m 

/(O, X, v) = (f ^) (x, v), (x, v)eM.^x (10) 
where the averaged relaxation operator is given by 

{QB)fix,v) = ^ [ [ Si\vlvs,\'^lv',,z){Miv)f(^,xs,v')-M{v')fix,v)}dv'dx[dx', 
with z = UcX + -^v — {oJcX' + -"-f ') and the averaged scattering cross section writes 
S{r, V3, r', W3, z) = o-(a/|z|2 + {v^ - ) x(r, r', z) 



with 



l{|r-r'|<|2|<r+r'} ^ ^ to, „, „,/ ^ to ^ ^ m2 



tt^a/I^^P — (r — r'Y's/ {t + ^')^ ~" kl^ 



The averaging technique allows us to treat many different collision operators, for ex- 
ample the Fokker-Planck kernel (see Appendix |A] for details) 

QFp{f) = — div, hj + ^vf) = — div, (mV. ({-) ] . 

mr V / niT { \M J ) 
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Theorem 1.2 The limit model when £ \ of (j3]), ([2]) with Q = Qfp is given by 

dj + ■ V^/ + ^3^x3/ + - (Es) d,J = {Qfp) /, (t, x, v)eR+xR^x (11) 

i3 m 

/(O, X, f ) = (f ^) (x, t;), (x, t;) G X (12) 
where the averaged Fokker-Planck operator and the diffusion matrix C write 

{Qfp) fix,v) = — div^^^r,^ \mC\/^^x,v ( 

mr V M 



E h 
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Notice that the averaged Fokker-Planck operator contains no derivatives with respect 
to X3 since the diffusion matrix C has only zero entries on the third line and column; 
averaging the Fokker-Planck operator leads to diffusion in velocity but also with respect 
to the perpendicular position coordinates. 

Our main motivation concerns the bilinear Fokker-Planck-Landau equation, more 
exactly how to average kernels like 

Qfpl(/, f){v) = div, ( / a{\v - v'\)S{v - v')[f{v')VJ{v) - f {v)V f {v')] dv' 

where a denotes the scattering cross section and S{w) = I — is the orthogonal 
projection on the plane of normal w, cf. |20J. Recall that Qfpl satisfies the mass, 
momentum and kinetic energy balances 

/ Qfpl(/, /)dt; = 0, / vQFPL{fJ)^v = Q, [ M!g^^^(/, /) dt; = 0. 

Jr3 J^3 Jr3 I 

Moreover it decreases the entropy flnf since, by standard computations, we obtain 

InfQFPLifJ) dv 

3 

1 f f aK^-^') A(V.ln/(^;)-V..ln/(t;'))P 



f{v)f{v') '' ' ^ dv'dv < 



2 ViMS \v - v\ 

We expect that the averaged Fokker-Planck-Landau operator satisfies the same prop- 
erties. Nevertheless we will see that all of them hold true only globally in velocity and 
perpendicular position coordinates. Indeed, the averaged collision kernel will account 
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for the interactions between Larmor circles (characterized by the center x+ '^v/oJc 
and the radius |l;|/|a;c|) rather than between particles. We show that the averaged 
Fokker-Planck-Landau kernel has the form 

(Qfpl) (/,/): = {QfplUJ)){^,v) 

= div^^^^^lul axf{x', X3, v')A+V^^^^^f {x,v) dv'dx[dx'2 




- div,^^^,^ <^ / / axfix,v)A Vu,,x',v'fix',X3,v')dv'dx[dx2 



(13) 



with z = UcX + -^v — {(jJcX' + """f'), cr = ^"(a/I-zP + (^^3 — ^3)^ )i X = xd^l? I"^'!) -2) and 

4 

crxA+(r, r\ W3, z) = ^ f (x, v, x', v') ® f (x, v, x', v') 



i=l 
4 



axA (r, fa, r', v'^, z) = ^ £tC{,x, v, x', v') (g) f (x', v', x, v) 

i=l 

for some vector fields (^*)i<i<4 and ei = 62 = —1,£3 =£4 = 1, see Proposition 15.91 
Actually have only zero entries on the third line and column and therefore, 

averaging the Fokker-Planck-Landau kernel leads to diffusion (and convolution) with 
respect to velocity but also perpendicular position coordinates. To the best of our 
knowledge, this is the first completely explicit result on this topic. In particular, the 
above coUisional kernel decreases the entropy / In / since, by standard computations 
we obtain (see Proposition IS.lOp 

/ f\nf{Q^pL){fJ)dvdx,dx2 = -^J2f II hf 

X (f ■ Vln/ - ei{C)' ■ Vlnff di;'dx;dx2 di;dxidx2 < 0, X3 G M. 

Here, for any ^ M^, the notations ^ ^ rj stands for the matrix whose entries are 
(^ ® ri)ki = ^hTji, 1 < k,l < 6. We obtain formally the following stability result 

Theorem 1.3 Let us consider f"" > 0, (1 + | Inf °|)f ° G L\R^ x M^) and denote by 
the solution 0/ (EJ, ([2]) with Q = Qfpl, for any e > 0. Then the limit f = lim^x^o /^ 
satisfies 

dtf+^^-V^f+V3d^J+^{E3)d,J = {QFPL){fJ), {t,x,v)eR+xR'xM.' (14) 
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/(O, X, v) = (r) (x, v), (x, v)eR^ X M' 
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(15) 



where the averaged Fokker-Planck-Landau operator is given by fll3p . 

Our paper is organized as follows. In Section H] we introduce the average opera- 
tor along a characteristic flow. Section 12] is devoted to the commutation properties 
between average and first order differential operators. The average of the linear Boltz- 
mann kernel is computed in Section |H We establish its main properties and we prove 
the convergence result stated in Theorem 11.11 Section E] is dedicated to the bilinear 
Fokker-Planck-Landau kernel, Theorem 11.31 We give a explicit form of its average and 
check the main physical properties. We prove the mass, momentum and total energy 
conservations for smooth solutions of the averaged Fokker-Planck-Landau equation 
coupled to the Poisson equation for the electric field. We also show that the mean 
Larmor circle center and power (with respect to the origin) are left invariant. Up to 
our knowledge this has not been reported yet. 

2 Average operator 

We recall briefly the definition and properties of the average operator corresponding 
to the transport operator T, whose definition in the L^(]R^ x M.^) setting is 

Tu = div^. „(m b), b = (y, 0, Uc ^v, 0), ujc= — 

m 

for any function u in the domain 

D(r) = {u{x,v) e l2(m3 X . div^,^(M b) e l2(m3 X m?)]. 

We denote by || ■ || the standard norm of L^(R3 x R3).The characteristics (X, V){s\ x, v) 
associated to v ■ + coc'^v ■ Vv, see ([7j), satisfy 








ds 







implying that 



V 



V{s) = R{-uJcs)v, X{s) = x + 



Xsi'S) = X3, 1/3(5) = V3 
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where R{a) stands for the rotation of angle a 

R{a) = 



cos a — sm a 
sin a cos a 



All the trajectories are Tc = Itx juj^ periodic and we introduce the average operator, see 
[2], for any function u E L2(M^ x M^) 



1 '•^^ 



{u){x,v) = 7^ u{X{s;x,v),V{s;x,v)) ds 

-'c Jo 

f- ^{R{a)v} \ 
/ u\x^ , X3, itfajf , f3 da. 



f.27r 

27r 



It is convenient to introduce the notation e**^ for the vector {cos (p, simp). Assume 
that the vector v writes v = \v\e'^'^. Then R{a)v = \v\e^^°'~^'^^ and the expression for (u) 



becomes 



{u){x,v) = — / ulx + — -^^ ^,X3,|l7|e'("+^),t;3 da 

27r Jq \ Uc uJc J 

1 /"^^ /_ ^{\v\e"^} 



M hrH ,X3, |t;|e*",t;3 da. (16) 

2vr Jq \ Uc uJc J 

Notice that {u) depends only on the invariants x + \v\,X3, V3 and therefore belongs 

to kerT. The following two results are justified in |3j, Propositions 2.1, 2.2. The first 

one states that averaging reduces to orthogonal projection onto the kernel of T. The 

second one concerns the invertibility of T on the subspace of zero average functions 

and establishes a Poincare inequality. 

Proposition 2.1 The average operator is linear continuous. Moreover it coincides 
with the orthogonal projection on the kernel ofT i.e., 

(u) G kerT and / (u — {u))(p dvdx = 0, V (/? G kerT. (17) 

JR3 J ^3 

Remark 2.1 Notice that {X, V) depends only on s and {x,v) and thus the variational 
characterization in (|T7|) holds true at any fixed (X3, ^3) G M^. Indeed, for any (f G ker T, 
(x3,f3) G we have 

{uLp){x,v) dvdx = — I II u{x,v)ip{X{—s;x,v),X3,V{—s]x,v),V3)dvdxds 




JO 

1 '■^^ 



Jo 





u{X{s; X, f ), X3, V{s; x, v),vs)(p{x, v) df dxds 
u) (x, v)ip{x^ v) dUdaJ. 
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We have the orthogonal decomposition of L^(]R'^ x M.^) into invariant functions along 
the characteristics ([7]) and zero average functions 



u = (u) + {u — (m)), / / (m — {u)) {u) dvdx = 0. 
Notice that T* = — T and thus the equality (■) = Projj^gj.^- implies 



ker (■) = (ker T)^ = (ker T*)^ = Range T. 

In particular Range T C ker (■). Actually we show that Range T is closed, which will 
give a solvability condition for Tu = w (cf. [3J, Propositions 2.2). 

Proposition 2.2 The restriction of T to ker (■) is one to one map onto ker(-). Its 
inverse belongs to £(ker (■) ,ker (■)) and we have the Poincare inequality 

II II ^ 27r qB 

\\U\\ < r / U\\, iOc = 7^ U 

\Uc\ m 

for any u G D(T) fl ker (■). 

The natural space when dealing with the linear Boltzmann kernel Qb is L"^ {M~^dxdv) 
rather than L^{dxdv). Motivated by that we introduce the operator Tm '■ D(Tm) C 
(M~-^da;df) — (M~^dxdf ) given by 7a/m = diYx,v{ub) for any function u in the 
domain 

D{Tm) = {u{x,v) e L"^ {M-^dxdv) : dw^^^{ub) e L"^ (M-^dxdv)}. 



Straightforward arguments show that u G D{Tm) iff u/y/M G D(T) and Tm{u) = 
y/MT{u/\/M) for any u G V{Tm)- In particular we have kerT/if = v^MkerT. Notice 
that formula f|T6|l still defines a linear bounded operator on L'^ {M~^dxdv), denoted 
by i-)]^, which coincides with the orthogonal projection on the kernel of Tm, with 
respect to the scalar product of {M~^dxdv). Indeed, taking into account that M{v) 
is constant along the characteristic flow of ([7]), we have for any u G (M-^dxdv) 

{u)^ = y/M(-^) eVm ker r = ker71^ 



\ VM/ 
and for any ip G ker 7m 

dxdv f f f u / u \\ (p 
[u - {u}j^) ip{x, v) — — = i —= - ( —= ) —= dvdx = 0. 




M{v) JM.JM.3\^/M 1 ) 4M 
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The Poincare inequality holds also true, with the same constant, since for any u G 
D(Ta/) n ker (■)^^ we can write 

||^|U2(A./-i) = 

From now on, for the sake of simplicity, we will use only the notations T, (■), indepen- 
dently of acting on L^(da;df ) or (M^^dxdf ). 



u 



< 



2-K 



\UJr 



r 



u 



27r 



27r 



WTmuW 



L2(M-i)- 



3 Average and first order diff^erential operators 

We intend to average transport operators, see ([H]). Moreover, in order to handle the 
Fokker-Planck-Landau kernel we will need to average second order differential opera- 
tors. For doing that it is convenient to identify derivations which leave invariant ker T. 
It turns out that these derivations are those along the invariants 

We introduce also ipo = with v = \v\e^'^, a G [0,27r[. Notice that ipo has a jump 
of ^ across v G R+ x {0} but not its gradient with respect to v 

VvOt = -TUpj, ^vll>0 = PTJ, TlpO = 1. 

The idea is to consider the fields (&*)o<i<5 such that 

Indeed, the map {x,v) — )■ {ipi{x,v))o<i<5 defines a change of coordinates 
xi = i)i + — sin(cUcV'o), X2 = ^jj2 + — cos(wcV'o), X3 = V'3 

= cos{ujcipo), V2 = -ipi sin(u;c^o), ^^3 = V^5- 

Therefore any function u = u{x,v) can be written u{x,v) = U{iIj{x,v)), ip = {ipi)o<i<5 
and thus, for any i G {0, 1, 5} we have 

^ dU dU 
b' ■ Vx,vU = h'-Y^ -—{i){x,v))V^,^ipj = —{'4j{x,v)). 

j=0 

In other words the derivations ■'Vx,v act like d^^, < i < 5. In particular if m G ker T, 
meaning that U does not depend on ipo, then b^-'Vx,vU = d^.U^ip^x, v)) does not depend 
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on ipQ, saying tliat kerT is left invariant by 6* ■ Vx,v, < i < 5. The following result 
comes by direct computation and is left to the reader. For any smooth vector fields 
^,77 on MP, the notation [^,r)] stands for their Poisson bracket i.e., 



Proposition 3.1 The fields (&*)o<i<5 satisfying U ■ Vx,v'4^j — 5], i, j ^ are given 
by 



V3- 



\v\ 

Moreover the Poisson brackets between (&*)o<i<5 vanishes or equivalently the derivations 
b^ • ^x,v) < i < 5 are commuting. 

Remcirk 3.1 Notice that {b^)i^4 are divergence free and divx,vb'^ — 

We claim that the operators u diVa;,„(M6*), with domain 

D(div,;,^(- b')) = {ue L\M.^ X M^) : div^,^(M6^) G L\M.^ x R^)}, < i < 5 

are commuting with the average operator. More generally we establish the following 
result. 

Proposition 3.2 Assume that the field c ■ Vx,v is in involution with b ■ Vx,v — v ■ + 
uJc '^v • i.e., [c, 6] = 0. Then the operator divj;^^(- c) is commuting with the average 
operator associated to the flow ofb- Vx,v that is, for any function u e D(diVa;,^(- c)) its 
average {u) belongs to D(diVa;,^(- c)) and 

diyx,v{{u) c) = (div^_„(Mc)) . 
Proof. Let us consider u e D(diVj;,^(- c)). For any (p e C^(R^ x R^) n kerT we have 
/ / {diVx,v{uc)) dvdx — / / diVx,v{uc)cp dvdx — — / / uc ■Vx,vV dvdx. (18) 
But T(c • Vx,v^) = c ■ Vx,viT^) = saying that c ■ Vx,v^ G ker T and thus 

/ uc- Vx,v'^ dvdx —I I (u) c- Vx,v'^ dvdx. (19) 
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Combining ([18]), (IT9]) we obtain for any ^ G CliM? x B?) fl ker T 

(diVa;^„(Mc)) dwdx = — (u) c- Vx,v^ dvdx. (20) 




Actually the previous equality holds also true for smooth functions G ker (•). Indeed, 
by Proposition 12. 2[ for any smooth function (p G ker (■) there is ip & D(T) flker (•) such 
that TV" = and thus c ■ Vx,v'P = c ■ Vx,v(Tip) = T{c ■ Vx,vi') £ Range T = ker (■). 
Using now the orthogonality between kerT and ker (■) we deduce that 

/ / (div^,^(Mc))v9 dt;dx = = - / / {u) c-Vx,vV> dvdx, (/^ e C](]R^ x M^) nker (■) . 

Finally fl20p is verified for any smooth (p, implying that 

{u) G D(div^^^,(- c)) and divx,v{{u)c) = {divx,v{uc)) . 

□ 

We want to average transport operators, which are written in conservative forms. In 
order to obtain averaged model still written in conservative form, it is worth to establish 
the following commutation formula between average and divergence. For the sake of 
simphcity we discard all difficulties related to the required minimal smoothness. 

Proposition 3.3 For any smooth field ^ = {Cx,Cv) ^ we have the equality 
(div^,^^) = div5;|(^3r+— ■ ^ ^ 



U}(, / \ \V\ I i>Jr\V\ \ If I / UAV 







V 









+ div^l^e^-^^^ + ^e^- — ) — 1^ + 9,3(^.3) 

In 'particular we have for any smooth field C,x ^ 

(div^^^.) = div^ {^x) 

and for any smooth field G M.^ 

(div.^^) = divg. 1 ( — ) + (^v ■ -7^) —7=7 - i^v 

{ \ Uc / \ \V\ / Uc\v\ \ \V\ I UJc\V 
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Proof. By construction we have ^^=o ^* ® x,v'4'i = I and thus 

5 

The main statement follows thanks to Proposition 13. 2[ since we have 

(diva;,„0 = ■ ^x,v'^PiW]^ = div^^^ 1^ ■ Vx,v'ipi) . 

The other statements come by considering the fields (^x,0) and (O,,^^,). |^ 

A direct consequence of Proposition 13.31 is the computation of the average for the 
transport operator in ([6]). 

Proposition 3.4 Assume that the electric field derives from a smooth potential i.e., 
E = Then for any f G Cl(R^ x M^) p kerT we have 

(dtf + vsd^J + ■ V./ + Tf') = dtf + ■ V^/ + vsd^J + ^ {Es) d^J. 

\ m / B m 

Proof. We can write 

(dtf + v,d,J + ^E ■ V,/ + Tf') = dtf + {v.d^J) + ^{E- V,/) 

since {T f^) = and {dtf) = dt (/) = dtf . The average of v^d^^f comes easily thanks 
to Proposition 13.21 

Observe that T{f(p) = f v ■ = — / v ■ E and thus (^f v ■ E'^ = 0. Thanks to 
Proposition 13.31 one gets 

{divM'E}) = div^//— \ + r//^^\ + a,3 ifE,) 

\ LOr- / \ ujjvr / 



implying that 

^ (div,{/^}) = div^ I 1 + im- 

m \ B \ m 

Using again Proposition 13.21 notice that 

(^3) = div,,4(E3) 6^} = (div,,4E365}) = (9,3^3) = 

and 

d:iY^{^E) = {dw^^E) = 

and our statement follows. 
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□ 



Remark 3.2 We have proved that averaging the transport operator a ■ Vx,v '■= v^d^g + 
■ Vy leads to A ■ V^^y := ^-^-^ ■ Vgr + v^d^^ + ^ {E3) dy^ which verifies 

{a ■ V.,./) = A ■ V.,yf, f G Cl {R' X M^) n ker T. 

By construction, the operator A ■ Vx,v leaves invariant the suhspace of smooth functions 
of ker T ■ By antisymmetry (since div^^^A = 0) it is easily seen that A ■ Vx,y also leaves 
invariant the suhspace of smooth functions in ker (■). Indeed, consider h a zero average 
smooth function and let us prove that {A ■ Vx,vh) = For any smooth f in kerT we 
have 

/ A- Vx,yh f dvdx = - HA - Vx,yf dvdx = 

by the orthogonality between ker (■) andkeiT, and thus {A ■ Vx,vh) = 0. Finally A-'Vx,v 
is commuting with the average operator {A ■ Vx,vf) = A ■ Vx,v (/) for any smooth f . 

4 The relaxation collision operator 

In this section we analyze the hnear Boltzmann coUision kernel [23l [22] 

QB{f){x,v) = - [ s{v,v'){M{v)f{x,v') - M{v')f{x,v)} dv' (21) 

T Jr3 

where the scattering cross section satisfies 

s{v,v') = s{v\v), < So < s(f,f') < 5*0 < +00, f,f'GR^ (22) 

We recall the standard properties of this operator. Here denote the gain/loss 
relaxation collision operators 

QbU){v) = -I siv,v')M{v)fiv')dv\ Q-{f){v) = -[ s{v,v')M{v')fiv)dv'. 

Proposition 4.1 Assume that the scattering cross section satisfies fl22|) . Then 

1. The gain/loss collision operators are linear bounded operators on L'^{M~^dv), 
with WQ^W < Sq/t, and symmetric with respect to the scalar product of L'^(M~^dv) . 

2. For any f G L'^{M^^dv) we have 

I QB{f){v)f{v)^ = -^ f f s{v,v')M{v)M{v') 



f{v) f{v') 
M{v) M{v') 



dv' dv < 0. 
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We want to average Qsif) for functions / satisfying the constraint (|5]). In this section 
the operators T and (■) should be understood in the (M~^dxdf ) framework. We 
need to compute the average of functions hke f^^ C{v,v')f{x,v') dv' where C{v,v') is 
a given function. The corresponding result in the bidimensional framework has been 
announced in [6j. We will see that we only need to consider functions C which are 
left invariant by any rotation around 63 = (0,0, 1). Therefore we assume that for any 
orthogonal matrix O G A^3(]R) such that Oe^ = 63 we have 

CCOv, 'Ov') = C{v, v'), V, v' e Ml (23) 

These functions are precisely those depending only on t>3, |t>'|, and the angle 
between v and v' 

C{v, v') = C{\v\,Vs, \v'\, fg, ip), (p = argt>' — axgv. 



Proposition 4.2 Assume that the function C{v,v') satisfies f l23p and belongs to the 
space L'^{M~^{v)M{v')dvdv') . Then for any function f G kerT we have 

C{v,v')f{x,v')dv'\{x,v) = uj'^ / / C{\v\,V3,\v'\,v';^, z)f{x',X3,v') dv'dx'idx2 

(24) 

where z = UcX + -^v — {ucX' + ^v') 



and for any \z\ G {\r — r'\,r + r') , ip G (0, vr) is the unique angle such that 



|2 2 I / /\2 o ' 

z\ =r + [r ) — 2rr cos (f. 



Proof. For (x, f ) G x we have 



C{v,v')f{x,v') dv' ] < / {C{v,v')yM{v')dv 



M{v') 



dv' 



implying that 



C(f , v')f{x, v') dv' 



< \\C\ 



L2{M~^dxdv) 



L'2{M-'^{v)M{v')dvdv')\\J \\L^{M-'^dxdv) 
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Therefore the function {x,v) — )■ J^^ C{v,v')f{x,v') dv' belongs to {M'^dxdv) and 
can be averaged in (M~-'^da;df) . Consider {x,v) G M'^ x M^. By formula flTB]) we 
have 



I: = (^j^C{v,v')f{xy) dv'^ {x,v) 



= ^fl C{W^,v,,v')f(x+ — -^-^^^^,x^,v'] dv'da. 
For any fixed a G [0, 2n) we use the cylindrical coordinates 

v' = (rV('^+"),t;^), r' G M+, [-7r,7r) 

and therefore 

X f(x + — - tM^^ r'e*('^+"), v'^ r'dr'd^dv'r^da. (25) 
But / G ker T and thus there is g such that 

f{x,v)=g(x + —,X3,\v\,V3^ (26) 



c 



implying that 

By one hand notice that r'e*'-'''^"-' — |iJ|e*° = Ze**^'^'*'"^ where = + (r')^ — 2rr' cos^p, 
r = \v\ and ip depends on r, r', but not on a. By the other hand, since C is invariant 
by rotation around 63 we deduce that 

C(re*", f3, r'e*^'^"'""-*, ^3) = (7(r, fs, r', </)), = arg f ' — arg tJ. 



The map (p — > /(</)) = \/r'^ ^ ir'Y' — 2rr' cosip defines a coordinate change between 
ip G (0, vr) and / G (|r — r'|, r + r') and 

2/d/ 

dip 
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By Fubini theorem one gets 



I = — / C{r,vs,r',v'^,ip)g[x + — + ,X3,r' ,v'^] r' dadr'dipdv'^ 

TT / C{r,V3,r',Vs,(p)g[x+ — + ,X3,r',v's) r' dadr'dipdv's 

Jk+ Jo \ Uc oOc J 

1 




-TT JR+ Jo 

{C(r, vs, r', v'^, ip) + C{r, v^, r', v'^, -(p)} 

s+ Jo Jo 




27r 



X g [x -\ 1 , X3,r ] r d(pdadr dv^ 



1 

2^ 



2tt fr+r' 

{C{r, V3, r', v^, <p{l)) + C{r, v^, r', v'^, -^{l))} 

+ Jo J\r-r'\ 

V ^{/e'"} , A 2lr' dldadr'dv' 




xg\x^ \ ,X3,r 



3 



For any a' e [0, 27r) we have 



x H \ , X3,r ,v^] = / X H \ ,X3, r e'" , 

and performing the change of coordinates v' = (r'e*"', ^3) leads to 

-1 p2n p p p27v p 

^=V-2 / / / {^(r,^3,r',^^,^(/)) + C(r,^3,r',^^,-^(/))} 

Jo Jm.+ Jr Jo Jr+ 

J- , , ^{^e^°} ^{r'e'"'} , ,\ l{\r-r'\<i<r+r'}r'dr'da'dv'^ldlda 
X f [x^ \ 5^ -, X3, re*" , Vn , ' ' ^ 



27r2 




{C(r, W3, W3, +C(r, W3, |w'|,W3, -<^(0)} 

'0 -/m+ -/IR3 



^^c t^c a;c ' ' y _ (1^1 _ \y,\Y^^\y\ + |^/|)2 _ 
Using the notation 

where for any \z\ e (|r — r'|, r + r'), </? e (0, tt) is the unique angle such that 

\z\^ = + {r'Y - 2rr'cosip 

one gets 

"271 



f f f C{\vlv3,\v'lv's,-^{le'''})f (x + — + ^^^-^,Xs,v') Idldadv'. 
Jm3 Jo Jr+ \ J 
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We take as new coordinates 



x' = x-\ h — ^ 

Ur iOr. UJr 



Observing that det ^^an'^\^ = -r "we deduce that 




□ 



Remark 4.1 The constraint Tf = allows us to reduce the right hand side of 
to a four dimensional integral. Indeed, thanks to fl26p we obtain, using the notation 
y = x + ^ 

I = ujI / / C(|l;|, f3, |f 'I, fg, (wc^ + ^t;) — (wcX' + ^v'))g{x'-\ ^, X3, |f '|, tig ) df'dx 



I I C{\v\,v:i,\v'\,v'^,ujc{y-y')) g{y',x:i,\v'\,v'^)dv'dy[dy'2 
2ttujI I I I C{\v\,v:i,r\v'^,ujc{y-y')) g{y',x:i,r\v'^)r'dr'dv'^dy'. 



We prefer to keep the five dimensional integral representation since, in the sequel, we 
will introduce similar integral terms, hut with densities f not satisfying the constraint 
Tf = 0. 

Remark 4.2 // the function C{r, V3, r', Wg, (p) is odd with respect to (p, then C = and 

^ j C{v, v')f{x, v') dv'"^ = 0, / G ker T. 
Remark 4.3 Let x be the function 

'^{\r-r'\<\z\<r+r'} 



X{r,r ,z) 



tt^a/I-zP — (r — r')^^/ {r + r'Y — \z\ 



for any r,r' G M+,z G M^. Then for every r,r' G M+, ''^'^ -2)d;z is a probability 
measure on 

/ x(r, r', dz; = 1, r,r' ElS.^ 

and ^ Jjj3 C(f , f f ') df') appears as a convolution with respect to the invariants 
UcX + -^v. Indeed, using the formula 

— ^v^ .—. 



UJr. ■ ■ 3 



/(x',X3,t;') = g[x' + —,3:3, \v'\,v'^ 
19 



we obtain 

If C{v,v')f{x,v')dv'){x,v)= [ [ C{\^,V3,M,V'2,{UJ,X+ ^v)~{u,'^+ ^V)) 

X g [ x' -\ ,X3,|f'|,f3) d{ucx' + '^v')dv'. 



Remark 4.4 The conclusion of Proposition also holds true for hounded functions 
f which are constant along the flow ([7]), provided that C{v,v') G L°° {dv; {dv')) and 
satisfies ( 12^ . Indeed, in this case f f^^ C{v,v')f{x,v') dv' is bounded on L°°{dxdv) 



C{v,v')f{x,v') dv' 



< ||C'||L°°(dt;;Ll(dv')) II /IU°°{dxdi;) 

L°°{dxdv) 



and using the L°° version of the average operator, the same computations as those in 
the proof of Proposition \4-S\ show that 

(/ C{v,v')f{x,v')dv'){x,v)=ujl / C{\v\, V3,\v'\, v'^,z)f {x',X3,v') dv'dx[dx2. 

Corollary 4.1 Assume that the scattering cross section satisfies (12^ and 

s{v,v') = a{\v -v'\), v,v' eW^ (27) 

for some function a : IR+ — )■ ]R_|_ . Then for any f G ker T we have 

{Qsf)ix,v) = '^ [ [ Si\v\,vs,\'^lv',,z){Miv)f(^,xs,v')-Miv')fix,v)}dv'dx[dx', 

with z = UcX + -^v — {(jJcX' + -^v') and 



S{r,V3,r',v'3,z) = ^(a/I^;^ + {vs - v'^y )xir,r',z). 

Proof. Clearly the function C{v,v') = a{\v — v'\)M{v) satisfies (123 p . belongs to 
L'^{M-^{v)M{v')dvdv') and we have 

5(r, t>3, r', t>3, ip) = a{^/r'^ + (r')^ — 2rr'cosv9 + (fs — ) 

S{r,V3,r',v'^,z) = o-(a/|zP + {v^ - v'^)^ )xir,r',z). 
Thanks to Proposition 14.21 we obtain, with z = {ucX + -^v) — {ucX' + -^v') 

s{v,v')M{v)f{x,v')dv'j=u'^ / / S{\v\,V3,\v'\,v'^^ z)M{v)f{x',X3,v') dv'dx[dx2 
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Since / belongs to (M ^dxdv) and remains constant along the flow ([7]) we have 
s{v, v')M{v')f{x, v) dv') = fix, v)( [ s{v, v')M{v') dv' 



where the flrst average operator should be understood in the (M~^dxdf ) setting 
and the second one in the L°°{dxdv) setting. Remark 14.41 applied with C{v,v') = 
s{v,v')M{v') G L°°{dv; L^{dv')) and the constant function 1 G L°°{dxdv) yields 

( / s{v,v')M{v') dv') =ujI f [ S{\v\,V3, m,v'3, z)M{v') dv'dx'^dx'^. 

\Jr^ I JR?JR? 

Therefore we obtain 

s[v,v')M[v')f[x,v) df'\ = ull [ S{\v\,V3,\17\,v'3,z)M{v')f{x,v) dv'dx'^dx'^ 



and our statement follows immediately. 



□ 



We intend to extend the previous averaged collision operator to all densities /, not 
only those satisfying the constraint T/ = 0. Think that, when simulating numerically 
these models, the particle density may not satisfy exactly Tf = 0, and thus we need 
to construct such a extension. One possibility consists to appeal to the decomposition 
/ = (/) + (/ — (/)) and to neglect the fluctuation / — (/), leading to the operator 



f^QBf--={QB (/)) 

= ^ /" /" S{\v\,Vs,\'^\,v'„z){M{v) if) (^,xs,v') - M{v') if) {x,v)} dv'dx'.dx'. 



for any / G L"^ {M ^dxdv) . Clearly Qb coincides with (QbI) for any / G kerT. 
Notice that for any (x,f), {x'^,v'^) the function 

(x', f') — )■ iS(|w|, f3, |f 'I, ^3, WcX + ^V — {UcX' + "'"f'))M(t>) 



depends only on Ucx' + v', \v'\ and therefore, thanks to Remark 12.11 we obtain a 
simpler form 



QBf = — / S{\v\,v,,\v%v',,z){M{v)f{x',X3,v') - M{v') if) {x,v)} dv'dx'.dx',. 

Nevertheless notice that it is not possible to remove the average in the loss part of the 
previous operator. Since Qb and (■) are linear bounded operators on (M~^dxdf ) 
we deduce that Qb is linear bounded on (M^^dxdf) . The properties of Qb come 
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immediately from the properties of Qb, cf. Proposition 14.11 For example we have for 
any / G {M-^dxdv) 




f 

i ivi 



1 

27 




{Qb if)) dvdx 

i IVl 




Qb if) ^ dvdx 

3 IVl 




s{v,v')M{v)M{v') 
Another possible extension is given by 



{f){x,v) {f){x,v') 



1 2 



M{v) 



M{v') 



dv'dvdx < 0. 



(Qb)/:=— / / S{\%v^MA^z){M{v)f{x',x^,v') - M{v')f{x,v)} dv'dx\dx'^ 

(28) 

for any f E L"^ (M~^dxdf) , which is very similar to the operator in ( l2T|) . We keep 
this operator as extension for the operator in Corollary 14.11 The properties of (Qb) 
are summarized below 

Proposition 4.3 Assume that the scattering cross section satisfies fl22|) . (|27|) . Then 

1. The operator (Qb) is linear bounded on (M~-'^dxdf ) and symmetric with respect 
to the scalar product of (M^-'^dxdf ) . 

2. For any f E L"^ (M~^dxdf ) we have 




(Qb) if )jj dvdx = 

i IVI It 




S{\vlv,,\v%v'^,z)M{v)M{v') 



f{x,v) f{x',X3,v') 



dv dx^dxgdfdx < 0. 



M{v) M{v') 

Proof. 1. The boundedness of the loss part follows easily since it is the multiplication 
by the bounded function (see Remark IT3|) 



to: 



T 




S{\v\,V3, \v'\,v'^,z)M{v') dv'dx[dx2 



T 




S{\v\,V3, \v'\,v'3,-z')M{v') dz'dv' 



r 



For the gain part we use the inequalities 



< 





S 



S{\v\,V3, \v'\,v'3,z)f{x',X3,v') dv'dx[dx2 
i 

P{¥,X3,V') 




s 



dv'dx[dx'2 ul I [ SM{v') dv'dx\dx^ 
df'dx'^dx2. 



M{v') 

P{¥,X3,V ') j^,,^^,^^, 

M{v') 
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Thanks to Remark [4.31 we deduce, with = (M ^da;di;) , that 



1^: 




— / / SM{v)f{x',X3,v') dv'dx[dx'^ 



^, — ^ df'dx'idx2dvdx 




< ^ / / M(t;)S'o 




M(t;') 



5 / Miv) 




P{X',X3,V') 

M{v') 



i S da;ida;2dt> dX]^dx2dx3 dv 



< ^ 




2. Interchanging {x',v') with {x,v) and observing that this change leaves invariant S, 
yield for any f,g & (M~-'^dxdf ) 



i{QB)f,g)Li 

,2 



r 




r 




SM{v)M{v') 



SM{v)M{v') 



f{x',X3,v') _ 

M{v') M{v) 
f{x',X3,v') _ f{x,v) 
M{v') M{v) 



dv'dx\dx'^ ^[ ^] } dvdx 
^ ^ M{v) 

dfdxida;2 



X , , ^' — ^df 'dx'^dxgdxs 



2 

c 
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5M(f)M(t;') 



M{v') M{v) 



g{x',X3,v') _ g{x,v) 
M{v') M{v) 



dv dxidx2 dvdx 



which justifies the symmetry of {Qb) and its negativity. 



□ 



Remark 4.5 Contrary to Qb, the operator {Qb) is non local in space. The value of 
{Qb) f o,t the point {x,v) depends on all the values of f in the set 

A{x,v) = {{x[,x'2,X3,v') : S{\v\,V3,\v'\,v'2,{ujcX + -^v) - {UcX' + ^v') ) > 0} 

= {{x[, X2, X3,v') : \ \v\ — \v'\ I < \{UcX + ^v) — {UcX' + ^v')\ < \v\ + \v'\}. 

Observe that if we denote by C^^v the Larmor circle 

Cx,v = {{x[,X2,X3) : I UJ^X' - {ujcX+ ^v) \ = \v\} 



then we have 



C^^y X {v' : v' e R^} C A{x,v) 



where X stands for the adherence of X in R^. In particular {x} x C A{x,v) 
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Remark 4.6 The gain/loss parts of (Qb) ore bounded on L^(M^ x M^) and for any 
f G L^(]R^ xM.^) we have the global balance of the mass J^sJ^aiQB) f dvdx = 0. Indeed, 
we have 

\\{Qb)^ f\W< — I I I I S{\^,v^M,v'^,z)M{v)\f{^,x^,v')\dv'dx\dx'^dvdx 
<— [ [ [ M{v)\f(x',X3,v')\ dv'dx^dx'^dxsdv 

= ^\\f\\L^ 

and similarly 

II (Qb)-/||li<^||/||li. 

r 

The global balance follows by interchanging {x',v') with {x,v). 

Combining (jS]), Propositions I3.4[ 14. II and ( 12 8 p yields the limit model in Theorem ll.il 
Proof, (of Theorem 0]) Clearly < G (M+,L1(M3 x M^)) and 

/ f'{t,x,v) dvdx = / / f'''{x,v) dvdx, teR+,£>0. 

We consider a sequence {ek)k C M* converging to such that limjt_>.+oo f^'' = f weakly -k 
in L°° (M4., (M~^dxdf )) . Using the weak formulation of ([3]), ([2]) with test functions 
rj{t)ip{x,v), Tj G C^(M+),(y9 G C^(]R^ X M^), we deduce, after multiplication by Sk and 
letting k 00, that the limit density satisfies the constraint 

r/(t) = 0, tGM+. (29) 

Considering test functions like r]{t)(p{x,v) with r] G Cc^(]R+), ip G Cg^(]R^ x M^) fikerT 
one gets 

/ / f"'W<^ + VV3dccs'^ + r]—E{x)-V^'^}dvdxdt+ f'"'r]{0)^ dvdx 

Jr+Jr3Jm_3 m JrsJr3 

= - [ [ [ VQBif')^ dvdxdt. 
jR+jR'ijR^ 

(30) 
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The symmetry of Qb cf. Proposition 14.11 allows us to write 

lim / / [ r]QB{r')^dvdxdt= lim / / [ r^^QBivM) 

fc^ + OO /rp , /in,3/in.3 k^ + OO J m _ m3 J m3 IVl 



[ [ vfQsivM) ^dt = [ [ [ vQbU)v dvdxdt 

Jr3JR3 IVl JR+JR3JR3 

r, {QbU))v dvdxdt = f [ f ri{QB)U)v dvdxdt 




(31) 

since f{t) G kerT.t G M+ and thus {QbU)) = (Qb) (/)• For the other terms in ([30]) 
we obtain thanks to Proposition 13. 4[ Remark 13.21 

lim / / / f ^ {dt + a- V,,,){r]^) dvdxdt = [ [ [ f {dt + a ■ V,,,){r]yD) dvdxdt 

k^+oo J^_^_J^:iJ^3 Jr+Jr3Jr3' 

f{dt + A-V,,,){r]^) dvdxdt 




(32) 

and 

f °(x, v)r]{0)<^{x, v) dvdx =11 r/(0) (f °) (x, v)'^{x, v) dvdx. (33) 




Combining ( 15U]) . (122]), dSS]) yields for any smooth test function ri(t)ip{x,v) with 

(y9 G ker T. 

/ / [ f{dt + A ■V,,,){w) dvdxdt + [ [ {r)v{OMx,v)dvdx 

^R+JR3jR3 JR3Jr3 

= - [ [ f {Qb) U)w dvdxdt. (34) 

JR+ JR3 JlR3 

By Remark 13.21 we know that A ■ V^,?; leaves invariant the subspace of zero average 
functions and therefore it is easily seen that flMl) is trivially satisfied for any test 
function r](t)i'{x,v), with G C](]R^ x M^) fi ker (■). Finally (l34l) holds true for any 
smooth test function, saying that / solves ([9]), (|T0|) . We are done if we prove the 
uniqueness for the solution of (Q, (fTOj) (and in this case all the family {f'^)e will 
converge towards /, weakly -k in L°° (]R+,L^ (M^-'^dxdt')) ). Assume that / solves 
with zero initial condition. By standard arguments one gets 

^■V^|/|+^;39.3l/| + - 
B m 

implying that 
d 



5*1/1 + ■ V^/l + vsd.,\f\ + ^ (E,) I/I = (Qb) (/) sgn/ 



dt 



\f{t,x,v) \ dvdx = (Qb) (/) sgnf{t,x,v) dvdx, t G 

' JR3JM3 
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Our conclusion comes by observing that 
/ / {Qb) (/) sgn/ Mx 

S{M{v)f(t,x',X3,v') — M{v')f(t,x,v)}sgnf(t,x,v) dt>'da;^da;2df dx 

3 

SM{v'){f(t,x,v)sgnf(t,x',X3,v') — \ f(t,x,v)\} df'dx']^dx2 dfdx < 0. 



T 

, ,2 





□ 

Remark 4.7 It is easily seen that the integro- differential operator in ([9]) propagates 
the constraint Tf = 0. We are done if we prove that fs = f for any s G M, where 
fs(t,x,v) = f(t,X{s;x,v),V{s;x,v)) and {X,V) is the characteristic flow (jTj). A 
direct computation shows thatT and A-Vx,v = ^ "V^ + f 36)^3 + ^ (-E3) d^^ commute, 
implying that 

^ ' ^x,vfs = ■ ^x,vf)s- 

Observe also that 

(Qb)^ fs = {Qb)^ f = {{Qb)^ f)s, {Qb)- fs = {{Qb)- f)s 

and therefore {Qb) fs = {{Qb) f)s- Finally both f,fs satisfy (Q, (|T0|) and our state- 
ment follows by the uniqueness that we have established before. 

Clearly the transport equation (Q can be written in the reduced phase space {y = 
X + -^5X3, r = |w|,f3) since, by the constraint Tf = 0, we know that f{t,x,v) = 
g{t,y,X3,r,V3) . We obtain 

9t9 + ^^^-'^y9 + Vsdx,g + -{Es)d,,g = 27T^ [ [ [ S{r,vs,r',v',,u,{y -^)) 

X {M g{t, y', X3, r', v'^^) - M' g{t, y, X3, r, V3)} r'dr'dv'^dy' 

where 

M = i e-^('-'+(''^)') M' = i e-t(('-')'+K)') 

(27r^/m)3/2 ' (27re/m)3/2 

Remark 4.8 The family (/^)o<e<i remains bounded in (]R_|_, (M~^dxdf )) for 
potentials of the form (f){x) = + (^^{xs). Indeed, in this case observe that the 
energy function W^{x, v) := ^^^^^ + q{e (j){x) + (f)3{x3)) satisfies 

dtW' + -{v-V^ + ujc^v- Vtj)W' + v^d^AV' + —Eix) ■ V^W' = 
£ m 
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and therefore one gets 

1 a ^ f if^^"^ 



^ Qs^nr ' 

M(t;)exp(-|[£ + 03(0:3)]) ■ 

Integrating with respect to (x, v) yields the hound 

(^'('' ^- d„d. < / / tel^ d.„d. 




lw?J^3M{v) exp(-|[£: 0(x) + 03(x3)]) Jr3J^3M{v) exp(-|[£: 0(x) + 03(3:3)]) 

implying the uniform estimate 

Ufsu ^ f \(l\ iUh^ + \\h\\L^) \ II .inn n^^^l 
11/ ||L°°(R+,L2(M-idxdt))) S exp I I 11/ \\L^(AI-i-dxdv), U < £ S i- 



5 The Fokker-Planck-Landau operator 

In this section we focus on the Fokker-Planck-Landau equation [TOj [TT| [T2| [T3] . The 
rate of change of the density fg, corresponding to a population of charged particles of 
specie s, due to colhsions with charge particles of specie s' writes 



QFPLifs,fs') = —diVv I ^ils'(^ss'{\v - V'\)\V - v'\^ 



X Siv-v') (^f^,(v')iVJs){v) - —fs{v){V.>fs'){v')] dv' 

where fiss' = m^+m , reduced mass of the pair {ms,ms/}, ass' = (^s's > is 

the scattering cross section between species {s, s'} and the matrix S{w) = ^/ — ^j^^ 
corresponds to the orthogonal projection on the plane orthogonal to w. As the electron 
mass is much smaller than the ion mass, we consider only the collisions between ions, 
whose distribution function is denoted by and satisfies 

(35) 

where g > is the ion charge and m is the ion mass. As in the relaxation case, 
the limit model comes by averaging the collision kernel Qfpl- The treatment of the 
Fokker-Planck-Landau kernel is much elaborated. Therefore we content ourselves of 
formal computations. The main properties of the Fokker-Planck-Landau operator are 
summarized below 
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Proposition 5.1 Consider the Fokker- Planck- Landau kernel between ions 
QfplUJ) = 7div. / a,,{\v-v'\)\v-vfS{v-v') {f{v'){VJ){v) - f{v){V.>f)iv')) dv'. 



4 

Then the mass, momentum and kinetic energy balances hold true 
[ mQppLif, f) = 0, ! mvQFPLif, f) dv = 0, [ m^-^QppLif, f) dv = 0. 

JR3 Jr3 Jm.3 ^ 

Moreover the entropy production D := — J^^il + In f)QppL{f, /) dv is non negative 
D = l [ [ a{\v~v'\)\v-v'\f{v)f{v')\{v-v')A{VJnf{v)-V,^\nf{v'))\^ dv'dv>0. 



Proof. All statements come easily by integration by parts, observing that Aii{v, v') + 

Au{v\v) = Q,^\iemAu{v,v') = au{\v-v'\)\v-vfS{v-v'){f{v')VJ{v)-f{v^^^ 

and S{y — v'){v — v') = 0. ^ 

With the notation a{\v — v'\) = ^aa^v — v'\)\v — v'\^ the collision kernel becomes 
QpPLif, f)iv) = div, / ai\v - v'DSiv - v') {fiv')iVJ)iv) - /(i;)(V,./)K)) dv'. 



5.1 Preliminary computations 

The Fokker-Planck-Landau operator combines convolution and differential operators 
in V. Therefore its average can be determined by studying the commutation proper- 
ties between convolution and derivation with respect to the average. First we apply 
the commutation formula between divergence and average. Next we are looking for 
commutation between convolution and average. It is convenient to split Qfpl into its 
gain and loss parts Q$pp. We introduce the following notations, for any function g 
and vectors Wi , W2 

{9).s ■= (^^ 9{x, v'Mlv - v'\)S{v - v') dv'"^ 
(9, wi)as '■= i^j ^ 9{x, v')a{\v - v'\)S{v - v')wi dv'"^ 

{g,wi,W2)^s ■= i^j , 9{x^v')a{\v - v'\)S{v - v') ■.wi®W2 dv'^ . 
Let us establish some useful formulae based on Proposition 14. 2[ For any orthogonal 
matrix O G A^3(M) we consider the application {v,v') — ?■ S{'^Ov — ^Ov'). It is easily 
seen that 

S{'Ov - 'Ov') = 'OS{v - v')0. 
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Notice also that for any orthogonal matrix O G A^3(M) such that Oe^ = 63 we have 



{^Ov, 0) = 0) and {^^Ov, 0) = *0( 0) for any v E R\ 

Lemma 5.1 The following applications are left invariant by any rotation around 63, 
that is they satisfy 



S{v - v') : iy, 0) ® (y, 0), S{v - v') : (y, 0) ® ( ^v, 0), S{v - v') : { ^v, 0) ® ( ^v, 0) 
S : 0) ® (U, 0), ^ : (7, 0) ® ( ^U, 0), ^ : ( V, 0) ® (t;, 0), ^ : ( ^17, 0) ® ( ^I;, 0). 
Proof. For any rotation O around 63 we have 



SCOv - 'Ov') : i'Ov, 0) ® 0) = 'OS{v - v')0 : *C(tJ, 0) ® '0{v, 0) 

The other invariances follow similarly. |— I 
Therefore the formula in Proposition 14. 21 applies to all previous functions and we obtain 

Proposition 5.2 For any z E M."^ such that \r — r'\ < \z\ < r + r' we denote by 
(p G (0,27r) the angle satisfying + (r')^ — 2rr'cosip = Then for any function 

f G ker T we have, with the notation z = {ucX + -^v) — {ucX' + -^v') 




{f,{v,0),{v,0))^,=u^^ / / a( + (vs - v',Y )f{x', X3, v')x{\v\, \v'\,z) 



2 _ r (r-/cosy)- , 



(/, (v, 0), ( ^v, 0))^^ = (/, ( ^v, 0), (v, 0))^^ = 



3. 




{f,i^v,0),{^v,0))^=ul / / ai^\z\^ + (vs - v',y )f{x',xs,v')x{\v\, \v'\,z) 



><ir^-r4^^^P^Adv'dx[dx', 

\z\^+ (^3-^)2 I 
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{f,{v',0),{v,0))^^ = ul / / ai^\z\^ + ivs-v',y)f{x',Xs,v')xi\v\Av'\,z) 



f , rr'(r cos ip — r')(r — r' cos ip) , ,,,,,, 
X < rr cos 09 -— > dt> dxi dxn 

{ \z\^ + {vs - v'^y 

5. 

(/, (7, 0), ( ^v, 0))^^ = (/, ( ^7, 0), {v, 0))^^ = 



{f,{\',0),{^v,0))^^ = ul / a{^\z\^ + {vs - v',y )f{x', x,, v')x{\vl \v'lz) 



X { rr' COS LP — —-7, ; — )■ dv'dx\dx'n. 

Proof. We need to compute the functions C, C defined in Proposition 14.21 In each 

case we have 

1. 

[{v — v') ■ v]"^ 



C{v, v') = a{\v- v'\)S{v - v') : (y, 0) ® {v, 0) = a{\v - v'\) { \v 



\v — ti'P 



~ / r (r — j'j'^ COS u^'] 

C{r, V3, r', Wg, yj) = o-(^r2 + (r')^ - 2rr' cos + (f 3 - v'^y ) <r^ - , ,0 , . 

C{r, vs, r', t;^, z) = x(r, r', + (t^s - t;^)^ ) - ^^(^^3 -'^^)2 } • 

2. 

C(t;, t;') = (T(|t; - v'\)S{v - v') : (U, 0) ® ( ^U, 0) 

|f — v'\ 

~ I , ^ ( / I 12 I / ^ (r^-rr' cos (/?)rr' sin 

C(r,t;3,r ,t;3,(^) = -a(VkP + (^^3 - ^^3)' ) 1^|2 ^ _ ^/ )2 

C = 0. 

3. 



If — f'P 
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, , / I \ / rr~^ 7 TTo \ f 2 r^(r')^ sin^ 9? 1 

C(r,'U3,r ,^3,^) = xyr,r ,z)a[y/\z\^ + (^3 - v^)^ ) <!r - ^^^2 ^ (^^^ _ y ) ■ 



C{v, v') = a{\v - v'\)S{v - v') : {v', 0) (g) {v, 0) 

= a[\v — V \) <v' ■ V — — — ^ > 

{ \v — v'f ) 

~, / / s f / (rr' COS 09 — (r')^)(?"^ — rr' COS (zj) 1 
C(r, v„r,v„^)^a^rr cos ^ + (,3 _ ,.)2 1 

^, / / \ / / \ r / rr'(r cos 09 — r')(r — r' cos (ij) 
C(r, ^3, r , z) = r ,z)a< rr cos (/? , ,0 , . 



5. 



C(v, ^;') = - v'\)S{v - v') : (v', 0) ® ( ^tJ, 0) 
= 

^/ / / X / /rT5 7 7T?x f / • (rr' COS (p - (r'y)rr' sin ip\ 

C{r,V3,r ,V2,ip) = cT(VFr + {^3-^3) ) ^-rr sm(p 1^|2 ^ _ ^/)2 J 

C(r,V3,r',v's,z) = 0. 



C(w, w') = a{\v- v'\)S{v -v'):{ V, 0) ® ( 0) 

rr'cosifi - |^|2^^^^_^/-)2 | 
C(r, f3, r', ^3, 2;) = x(r, r', 2;)cr(A/|2;|2 + (^3 - u^)2 ) J rr' cos (p 



+ ("^3 - 1^3)^ 



□ 



We also need to compute the averages 

(/,(^,0)),5, </,(^U,0)>^^, </,K0)>^^, (/,(^<0)>^^. (36) 
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Notice that the functions aS{v — v'){v, 0), aS{v — v'){ -^v, 0), aS{v — v'){v', 0), aS{v — 
v'){ 0) writes 

for some scalar functions D, D' , D^, D'^ depending on fa, fg and (p, the angle 
between (y, 0), (f', 0). Performing the same steps as in the proof of Proposition 14.21 we 
obtain (see Appendix |B] for proof details) 

Proposition 5.3 Consider D = D{\v\,V3,\v'\,v'^,{p), D' = D'{\v\,V3,\v'\,v'2,(p) two 
functions and 

D{v, v') =Dv + D' 7, Dsiv, v') = D V3 + D' v'^. 
Then for any f G ker T we have 
1. 

( [ D{v,v')fix,v') dv'\ =ujI I [ V{\v\,V3, PU^, iLO,:x+ - {cu,F+ V)) 

X f{x',X3,v') df'dx'j^dxj 

where 

— Zl Z2 ^ 

uyr, V3, r , v^, z) = — —[D{r, V3, r', v'^, ^)re-'^ + D{r, V3, r', v^, -ip)re'^ 

+ D\r, V3, r', v'3, (^)rV(^-'^) + D'{r, V3, r', v'3, -(^)r'e-'('^-'^)] 
X X{r,r',z). 

2. 

( [ D3{v,v')f{x,v') dv')=ujl [ [ V3{\v\,V3,\'^\,v'3,{u,x+ ^v) - {u,'^+ ^17)) 

\ 7r3 / JR2 JR3 

X f{x',X3,v') dw'dx'^dx2 

where 

V3{r,V3,r',v'3,z) = ^[{D{r,V3,r',V3,^) + D{r,V3,r' ,v'3, -ip))v3 

+ {D'{r, V3, r', v'3, if) + £)'(r, V3, r', v'3, -'p))v'3]x{r\ r', z). 
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The angles (p,ip & (O? ^) such that = + (r')^ — 2rr'cosip, (r')^ = + + 
2r\z\ cosip. 

It is worth analyzing the case of even/odd coefficients D, D'. 

Proposition 5.4 With the same notations as in Proposition \5.3\ assume that the func- 
tions D,D' are even with respect to (f. Then we have 

1. 

V{r,V3,r',v'3,z) = [{D^ip) + D'{(p)) rcos^^ + \z\D' {if)]x{r,r' , z)j^ 

2. 

V3{r,V3,r',V3,z) = [{D{ip) + D'{ip)) V3 + D'{ip){v'3 - V3)]x{r,r' , z). 
Proof. 1. Clearly we have 

^[D{r,vs,r',v'^,ip) re''''' + D{r,V3,r' ,v'2, -<^) re'"^] = D{r,V3,r' ,v'^,ip) r(cos?/',0) 
and 

l[D'{r, V3, r', v'3, if) r'e'('^-'^) + D'{r, V3, r\ v'3, -ip) r' e-'^"^-"^^] = D' r'(cos(^/' - <^), 0). 

Consider now the triangle of vertices O = (0,0), A = {r,0),A' = r'e**^ in M^. The 
definitions for ifjtp assure that \z\ = \AA'\ and that ip is the supplement of the angle 
opposite to OA'. Applying the cosine theorem with respect to the angle opposite to 
OA one gets 

= [r'Y + |zp - 2r'\z\ cos(V^ - (f). (37) 
Combining with the definition of ip yields 

= 2\z\'^ — 2r'\z\ cos(?/' — Lp) + 2r\z\ cos-ip 

implying 

r cos ijj — r' cos{ip — ip) + \z\ = 0. (38) 

Finally one gets 

^z 

Vir, V3, r', v'3, z) = [D{(p) r cosip + D'{(p) r' cos{ip — ip)]x{r, r', z)—- 

^z 

= [{D{ip) + D'if)) r cosip + D'{ip){r' cos{ip — ip) — r cosip)]x{^: ^\ ^^Yl 
= [{D{ip) + D'{ip)) rcosip + \z\D'{(p)]x{r, r', z) — . 
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2. It follows immediately observing that 

1^3{r, V3, r', v'^, z) = {D{(p) V3 + D'{(p) v'3)xir, r', z) 

= [(D(v9) + ^'(v^)) V, + D'{^){v', - v,Mr- r\ z) 



□ 



Proposition 5.5 With the same notations as in Proposition \5.3\ assume that the func- 
tions D, D' are odd with respect to ip. Then we have 



1. 



~ ~ z 

V{r, V3, r', Wg, z) = -[D{ip) + D'{ip)] r sinip x{r, r', z) — 



z\ 



2. 

V3{r,V3,r',v'3,z) = 0. 

Proof. 1. We have 

^[D{r,V3,r',v'2,f) re"''^ + L'(r, W3, r>3, -y?) re'"^] = -{0, D{r,V3,r',V2,ip) rsinip) 
and 

^[D'{r,vs,r',v'^,ip) r'e''^^-^^ + D' {r^v^^r'.v'^, rV^^^-^)] = -(0,D' r'sin(V' - ^)). 

The sine theorem applied in the triangle of vertices O = (0,0),A = (r, 0),A' = r'e*"^ 
implies 

rsmif) = r' sin{ilj — ip). 

We deduce that 

~ ~ z 

V{r, t>3, r', t>3, z) = —[D{p) rsmip + D'{p) r' smijp — p)] xij, r', z)- 



\z\ 



\D{ip) + D\p))\ r sin-?/' x(r, r', z) —. 

Z/ 



2. Clearly we have b{p) + D{-p) = D'{p) + D'{-p) = and therefore V3 = 0. □ 
The averages in (l36i) come immediately appealing to Propositions 15.41 15 .51 



Corollary 5.1 With the notations in Proposition \ 5.3\ we have for any function f G 
ker T 
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5. 




|z| 

Proof. 1. We consider the function D{v,v') = a{\v — v'\)S{v — v'){v,0) = {D v + 

D' v', D'{v'^ — fa)) where 

~ , / r^-rr' cosy? 

D{r,V3,r ,v^,ip} = a [1 



y,2 _|_ — 2rr' cos + (f 3 — f 3) 



~ / / N -rr' cos 



+ {r'y — 2rr' cos + (173 — v'^y 
••Is 



2 / 

Thanks to Proposition 15.41 and the identity rcosip = _!_^n_£2^ obtain 



— rr' cos if \ , . , 



X rcosih + b --7^ ; — — dv dxidx. 



-^c / / (^fix[,X2,X3,v')x 

r4^^f^^^^^^f#^ df'dxidx', 
kl +(^^3-^^3) kl kl 



and also 



// b\v'^-V3)f{xy)dv'\=ujl f [ af{x[,x'„xs,v')x 

+ {v-i - fg) 



- rr' COS v? , , n j / j / j / 
^ U12 I (^3 - ^3) df dxidxg 



which justifies the first statement. 

2. We take D{v,v') = a{\v - v'\)S{v - v')(y ,Q) = (D v + D' 1/, D{v3 - v'^)) where 

(r')^ — rr' cosy? 



D{r,V3,r',v'^,(p) = a 



^2 _|_ — 2rr'cos<y9 + (^3 — ^3)^ 
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D'{r,V3,r',v'3,ip) = a - 
Notice that by dST]), (I38]) we have 

r cosip + \z\ = r' cos{iIj — (f) 

and in this case we obtain 

{Dv + b' V)f{x, v') dv') = ul 



{r') — rr' cos ip 



j,2 _|_ (^^/y _ 2rr' cos + {v^ — v'^y 



{r'y + 


z 


2 _ ^2 ^^/^2 _ 


rr' cos (f 


2 


z 




z 






Crfix[,X2,X3,v')x 



(r'Y — rr' cos 09 
rcosV- + 2 1 - I |„ . jy^ 



/ / o-/(a;'i,X2,X3,w')x 

3 

{r'f-rr'cos^ {y^-v'^f^z ,,,,,, 



df'dx'^dxg 




+ ('^3 — v'^)"^ \z\ \z\ 



dv dx;^dx2 



and 



^(^3 - v'^)f{x, v') dv') =ujI / Xg, 



(r'f - rr' cos ip , 

X -r-TT^ ^ 7^2 '-'^3 "~ "^S-^ "^"^ dX]^dX2 



justifying the second statement. 

3. We take D{v,v') = -a{\v - 0) = {D v + D' V^Dv^ + D'v'^) 

where 



D{r,V3,r',v'3,^) = -o- 



rr sm (p 



-D'. 



J.2 _|_ (jjy _ 2rr' cos + (f 3 — v'^Y 
By Proposition 15.51 we deduce that (^f^^ D{v,v')f{x,v') dv'') = 0. Therefore 



:/,(^^,o)> 



crS 



a{\v-v'\){^v,0)f{x,v') dv'). 



Applying now Proposition 15.41 with D = a, D' = we obtain 

/ / a-{\v — v'\)f{x,v')v dv'\ = uj"^ / / cr/x cos-?/'— df'dx']^dx2 

\ilR3 / jRiJ^i \Z\ 

and finally 

</,(^^,0)>^ 



'aS 




a fx r cos di^'dx'^dxg 

3 \z\ 




-rr' cos f a 'a 'a ' 
f M2 ^> ™ dx^dxg. 
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4. As before, by Proposition 15.51 we have 



and by Proposition 15.41 we obtain 



/ / a{\v — v'\)f{x,v')v' dv'\ = oo'^ / / a fx {r cos ifj + \z\)--^ dv'dx[dx'2 

\Jr3 I JK^jRa \z\ 



At the end one gets 

(/, ( ^7, 0))^, = -ul [ [ i^^'"'^ {z, 0) d^'da;;dx',. 

Jr^Jr3 \z\ 

□ 

The last average we will need is {f)^g = {J^^ f{x,v')a{\v — v'\)S{v — v') dv'Y By 
similar computations as those in the proofs of Propositions 14. 2[ 15.31 we obtain (see 
Appendix [Bl for details) 

Proposition 5.6 For any function f G kerT we have 




U).s = ^l / / <^W\A' + {vz-v',Y )f{x',x^,v')x{\vl \v%z)S{ {^z,v',-v,) ) dv'dx'.dx'^^ 



5.2 The averaged Fokker-Planck-Landau operator 

We are ready to determine the average of the Fokker-Planck-Landau kernel. For 
the sake of presentation we treat separately the gain and loss parts. Recall that the 
Fokker-Planck-Landau gain part appears as a velocity diffusion, where the diffusion 
matrix is a convolution in velocity 

Q^plU. f) = div. cx{\v - v'\)S{v - v')f{v') dv' V./(t;)| . 

The averaged Fokker-Planck-Landau kernel will keep the same structure, nevertheless 
diffusion and convolution have to be considered both in velocity and space perpendic- 
ular directions, as we have already observed in the relaxation case (see Remark 14.31) . 
The proof is postponed to Appendix O 
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Proposition 5.7 For any function f = f{x,v) satisfying the constraint Tf = we 
have 

X A+V^,.,./(x, v) dv'dx[dx'^} (39) 
with divujcx = ^^div^, V^,^ = 



{r'Y sin^ ^ (t.3 - v'sY ( {v, 0) ( ^v, 0) 



+ 



+ 



+ (t>3 — fg)^] \ \V\ 

r — r' cosip f (y,0) (-""U, 0)^ ^ f {^z,0) 

~\ i ' 



\v\ ' Ifl 



+ 



[r cos if — r) 



\v\ ' 



)(^3-^^3) A^^,0) {V,0)\ ((^3-^3)l^>-k|e3 



where z = {UcX + ^v) — {UcX' + -^v') and for any r, r' G M+, 2; G S'uc/i that \r — r'\ < 
\z\ < r + r' , the angle (f G (0, tt) is given by \z\'^ = + (r')^ — 2rr' cos v?- 

Remark 5.1 Clearly is symmetric and positive. Notice also that the vectors (63,0) 
and ( {z, 0), {—^z, — v'^) ) are orthogonal on 

(U,0) {^v,0)\ f{^v,0) {v,0)\ f{^z,0) \ {i^3-v'3)^,-\z\e3^ 



\v\ \v\ 
saying that 



,0,, 



A+{e,,0) = A+{ {z,0), {-^z,V3-v',) ) = 0. 
Actually we have for any z 

ker A+(r,W3,r',W3,2;) = span{(e3,0), ( {z,0), {-^z,V3 - v'^) )}. 

A similar result can be carried out for the loss part Qpp^ (see Appendix [C] for the 
proof). 

Proposition 5.8 For any function f = f{x,v) satisfying the constraint Tf = we 
have 



<Q]^PL(/,/)> = div.,.,.K / / a{y/\z\^ + {v,-v',y)f{x,v) x(M,kU) 

f{x[,X2,X3,v') df'dx^dxa} (40) 
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with 



rr' sin^ if - v'r^y fiv,0) i^v,0)\ fiv',0) 



+ (f3 — fg)^] \ \v 



I /I ' I /I 

\v \ \v \ 



r -r'cos<^ f iy,0) {^v,0)\ A^-^, 0) 



\z\ \ \v\ \v\ 
r cos if — r' f (t>', 0) ( 0) 



+ 



\Z\ \ \V'\ \V' 

rr' sin^ 09 / ( -^v, 0) (v, 0) 



-^,0) 



,0 



\v\ 



\v\ 



V,0) (t;',0) 



+ 



\z\^\z\^ + {v,-v',y \ \v 



+ 



\z\^\z\^ + {v,-v',y V 1^' 



=7^ + 



where z = {UcX + ^v) — {UcX' + -^v') and for any r, r' G M+, 2; G such that \r — r'\ < 
\z\ < r + r' , the angle (f G (0, tt) is given by \z\'^ = + (r')^ — 2rr' cos v?. 

For any x = (xi, X2), = {x[, x'2) G M^, t;, f ' G we introduce the fields 

r' sin (y9 (t>3 — fg) f (v,0) ( ^tJ, 0) 







r — r' cos V9 / (U, 0) (^t;, 0)\ f {^z,0) 



\v\ ' It;! 



e{x,v,^,v') = {axV/'- 

e{x,v,l^,v') = {axV^' 
e(x,v,¥,v') = {axV/^- 

^'^{x,v,x',v') _ (r'cos^-r)(f3-f^) /(%,0) (t;,0)\ ^ ((^3 - fDl^- -kks) 



r' sirup /( -'"W, 0) (f , 0) 



VkP + (^3-t;^)^ 



where r = \v\,r' = \v'\,z = {oOcX + '^v) — {oocx' + '^v'),a = a^\z\^ + (t>3 — fg)^, x = 
x{r, r', z) and ip G (0, tt) is given by |zp = + (r')^ — 2rr' cos y?, |r — r'| < |2;| < r + r'. 
Thanks to Propositions 15.7^ 15.81 we obtain the representation formula 

Proposition 5.9 Consider a function f = f{x,v) satisfying the constraint Tf = 0. 
Then 
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1. The averaged Fokker-Planck-Landau operator writes 



OO;' {QFPL{f,f)){x,v) = 

diyuj,x,v { y2 f^^'^ ^3, v')i\x, V, x', v') O f (x, v, x', v')Vu;,x,vf{x, v) dv'dx^dx'^ 

-div^,:,.,^, W / 52 f{x, v)C{x, V, x', v') O e^f (x', v', x, t;) V^^^'y/(x', xg, v') d^'dx^d^ 

(41) 

where ei = 62 = —l,es = = 1. 

2. The following properties hold true for any fixed X3 G M 

/ / {QpFiif, /)) {x, v) df dxidxa = 0, / v {QppLif, /)) {x, v) dt;dxidx2 = 

I I ^-^{QpFhifJ)) {x,v) dfdxidx2 = 0, /" I In f {QfplU'J')) i^^'") dfdxidxa < 
Proof. 1. By Proposition ISTTl we know that 

u-'{Q^FpdfJl)i^,v)= (42) 

div,^^^,^ W / 52 ^3' V, x', v')V^^x^^f{x, v) df 'dx'idx2 } . 

[Jr2Jrs~1 J 

Observe that we have x(r', r, — 2;) = x(r, r',z). Therefore the permutation {x,v) < — t- 
(x', v') leads to 

klvkl +(^^3-^3)^ V 1^^ I 1^^ I 



^^{x',v',x,v) = {ax} 



1/2 



r' -rcosip f{v',0) {^v',0)\ f{^z,0) 



\Z\ 



e'(x', t;', X, t;) = {axV^^^Y^ . - 

e^(^,^;^x,^) _ {rcos'f-r'){v's-v^) / ( ^7,0) (7,0) \ ((^3 - ^3)l|^> -kks) 

W}'/' ~ |;2|v/kP + (ti3-^)2 V H ' H V|z|2 + (t;3-t;^)2 

where z = {ucX + '^v) — {ujcX' + -^v'). By Proposition 15.81 one gets 

ujf{QFpdfJ)){x,v)= (43) 

div^^^,^ 1^2 ^)^*(^' ^' ® ^^^(a^', X, v)Vuj,x'.v'fix', X3, t;') dt;'dx^dx2 
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and the first statement follows combining (H2l) and fH3|) . 

2. The mass, third momentum component and kinetic energy balances for the averaged 
Fokker-Planck-Landau operator come by the corresponding properties of the Fokker- 
Planck-Landau kernel. Indeed, since l,f3,^ belong to kerT, we can write for any 
X3 G M, thanks to Remark 12.11 




1, ^3, -TT- (■ {QppLif, /)) dt;dxidx2 = / <l,V3, — > QppLif, f) dvdxidx2 

= (0,0,0). 

Similarly since In / G ker T one gets 

/ In f {QfplUJ)) dvdxidx2= \n f QfplU, f) dvdxidx2 (44) 
' ' ' a{\v — V \)j[x,vjj[x,v ) ^-j — df df dxidx2 



2 ilRZ JksJks F — "f , 

Finally observe that {vi) = {V2) = 0, {QppLif, f)) G kerT and thus trivially 

[vi,V2) {QfplU,!)) dt;da;ida;2 = (0,0). 




□ 



We establish formally the limit model stated in Theorem 11.31 

Proof, (of Theorem OJ Plugging the Ansatz = f + ef^ + e'^p + ... into (ESD we 
obtain 

(dt + vA, + -E-V, + -r] if + + ...) = QppLif, f) 
\ m e J 

+ e{QpPL{f,f)+QFPL{f\f)) + ... 



implying that 



r/ = 0, dtf + v^d^j + ■ V./ + rp = QppLif, /). 

m 



Applying the average operator, we deduce by Propositions 13. 4[ 15.91 that / satisfies 



dj + ■ Vsr/ + vsd^j + ^ {E,) d^J = {Qfpl) (/, /). 

B m 



□ 
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As for the relaxation Boltzmann operator, we are searching for extensions of the 
averaged Fokker-Planck-Landau operator to the whole space of densities / = f{x,v), 
not necessarily in the kernel of T. One possibility is to consider the extension (Qfpl) 
obtained thanks to f HTj) . that is for any / 

oo:' {Qfpl) (/, /) = 

div<^^^,„ W / E .^(^^ ^3, v')C{x, V, x', v') ® V, x', v')V^^^^^f{x, v) dv'dx[dx2 

-div^^^^^ W / 52 ^)^'(^' ^' ^0 ® v', X, v)V:^^^'yf{x', xs, v') dv'dx[dx2 

(45) 

What is remarkable is that this extension still satisfies the mass, third momentum 
component, kinetic energy balances and decreases the entropy /In/, globally in (x,v). 

Proposition 5.10 Consider two functions f = f{x,v),ip = ip{x,v). For any X3 G M 
we have 

{Qfpl) (/, f)^ dt^dxidxa = -^x (46) 




4 



i=l 

where 



H / , / / ff'^^' ■ Vln/ - Si{C)'V\nf){C ■ - SiiCyV^') dv'dx\dx'^ dvdx^dx^ 



/ = f{x,v), f = /'(x'i,X2,X3,w'), Vif = V^^^,MX,V), Vif' = V^^a:',v'^{x[,X2,Xs,v') 

In particular the averaged Fokker-Planck-Landau operator satisfies the mass, third 
momentum component, kinetic energy balances (globally in {xi,X2,v)) 

[ [ f 1, V,, ^) {Qfpl) (/, /) di;dxidx2 = (0, 0, 0) 
7r2 7r3V ^ / 

and decreases the entropy / In/ (globally in {xi,X2,v)) 

I / In / {Qfpl) (/, /) dvdx^dx^ < 0. 

7r2 Jr3 
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Proof. Notice that for any 1 < z < 4 we have C,^ ■ (es, 0) = and therefore the operator 
d^^ujcx,v acts only in (xi, X2,v). Thus, for any fixed X3 G M we can perform integration 
by parts with respect to {xi,X2,v). 

f [{QFPL){fJ)vdvdx,dx2 = ~f2''cf f f Iff (47) 
X {(f ■ V<^)(f ■ Vln/) - ei{C ■ V<^)((r)' ■ V'ln/')} dt;'da;;dx^ Mxi<\x2. 
Performing the change of variables {x'^^^x^iv') o (xi,X2,f) yields 

/ f{QFPL)ifJWdvdx^dx2 = -Y,^l[ f f Iff (48) 
X {((e)' ■ VV)((e)' • V'ln/') - £,((f )' ■ V'vp')(e ■ Vln/)} dvdxAx2 dt;'da;;d4. 
Combining fj47|) . f l48|) one gets by Fubini theorem 

/ / (Qfpl) if, f)f df dxidx2 = / / / / //T* dt;'dxidx'2 df dxidx2 

where 

r = (e-V^-£,(a'-V'^') (e-Vln/-£,(e)'-V'ln/'), l<z<4. 

Clearly, the divergence form of {Qfpl) guarantees the mass conservation and 
applied with (f = \n f ensures that the entropy / In / is decreasing 

/ [\nf{QFPL)if,f)dvdx,dx2 = -^J2[ II Iff 



X 



(f • Vln / - EiiCY ■ Vlnf'Y dt;'dx;dx^ dfdxidxa < 0, X3 € 



It remains to show the kinetic energy and third momentum component balances. 
Thanks to formula fH6|) it is sufficient to check that 

e-V^-£.(e)'-V'^ = 0, l<z<4. 

The above condition is trivially satisfied for i G {1,2}. For i = 3 we have 
*q ^I^^P ,^-^M „,k''P r 1 1 /9 '^' sin 09 . - -I /g r sin (z; , 
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Finally, when i = 4 we obtain 



2 ' 2 



11/2! {r' cosip — r){v3 — v'^)r + \z\'^V3 {r coscp — r'){v!^ — V'^)r' + \z\'^v'^ 
= WxV^' „ 7 [r^ + (r')^ - 2rr'cos^] = 0. 

vfI + (^3 -^^3) 

Notice also that for any i G {1, 2, 3, 4} we have 
saying that {Qfpl) satisfies 

^3 {Qfpl) (/, /) df dxidx2 = 0. 
Jr2 7r3 



□ 



Remark 5.2 By the formula (HBI) we deduce that the positive smooth functions f sat- 
isfying {Qfpl) (/, /) = are those verifying 

r ■ V In / - e.iCy ■ V In /' = 0, ^ e {1, 2, 3, 4}. (49) 

In particular fl49p /ioWs trwe /or any Maxwellian f which belongs to ker T, since in that 
case 

{Qfpl) (/,/) = (Qfpl(/, /)) = (0) = 0. 

We deduce that 

r ■ v<^ - e.iey ■ vy = o, ^ g {i, 2, 3, 4} (50) 

for any function ip{x,v) = a{x)\v\'^ + /3{x)-v+'j{x) satisfying Tip = 0, and in particular 
for the functions 



Xi H , X2 , X3, V3, \x\ + 2x ■ — 

U!r. OJr. UJr 



X H 



\V 



|2 
2~' 



Notice that x + is the center of the circle obtained by projecting the Larmor circle 
Cx,v onto {x[, x'2) and \x\'^ + 2x ■ ^ is the power of the origin (0, 0) with respect to the 
same circle. We obtain, thanks to (1461) and f l50l) . the balances 




x + —] {Qfpl) if, f) dt;dxidx2 = (0, 0) 
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|xp + 2x ■ — ) {Qfpl) if, f) dt;dxidx2 = 



for any smooth function f . 

The previous identities allow us to establish the mass, momentum, total energy, Larmor 
circle center and power conservations for smooth solutions of f|T^ . ( IT5|) coupled with 
the Poisson equation for the electric field 

E = -V^y?, Eo div^E(t, x) = q fit, x, v) dv =: pit, x), (t, x) G M+ x (51) 



Theorem 5.1 Assume that {f,E) is a smooth solution for 



dtf+^^-V^f+V3d^J+^ {E;)d,J = {Qfpl) (fj), Tf = 0, (t,x,t;) G M+xM^x^^ 
B m 



eodiv^E(t, x) = q fit, X, v) dv, E = {t, x) G M+ x 

Then we have 

— I I f dvdx = ~r [ [ dvdx = 0, — /" [ vf dvdx = 

d f^O /" I rn|2 1 , /" /" I^P rill « ^ [ [ /^I^P , — 



-<^^/ \E\'dx+ m^f dvdx \ =0, - [^ + x- — \fdvdx = 0. 

at I z j^3 Jr3Jk3 / J at j]K3j]R3\ 2 coc / 

Proof. The mass conservation comes by the conservative formulation of the Vlasov 

equation, which writes 



-E) 



dtf + div^ <^ ^ + d,,{fv,} + ^d,,{f {E,)} = {Qfpl) if, /)• (52) 

B \ m 

Multiplying (|52|) by v and integrating with respect to (x, v) yield 

^ ^ vf dvdx — f f —f {E'^) 63 df dx = 




dt 

since for functions in ker T we can write 



/ / V {Qfpl) if, f) dvdx = [ [ v{QFPLif,f)) dvdx = [ [ {v) Qppdf, f) dvdx 

JR3Jr3 Jr3JIR3 Jm.3Jr3 

iO,0, V3)QFPLif J) dt;dx = 0. 
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Using the Poisson equation and the identity 



div.££ = div.(£<S£)-iv.|£P 



(53) 



we obtain, taking into account that / G ker T 

[ [ —f {E-i) dvdx = — [ [ fEs dvdx 
i]R3j]R3m m J JUS J ^3 



£o 



E^ div^^-E" dx 



' ^diY,{EsE)-^d^,.jE\-'} dx = 



m 



and thus ^ /ms/ks"^/ dvdx = 0. Muhiplying (152|) by ^^^^y^ an^l integrating with respect 
to yield 



dt 



mm 



-/ df dx — / / qfvs (-E3) dfdx = 0. 



Thanks to the continuity equation 



dt I fdv + div^ / /i-^ dv + d,, I fvs dv = 



we have 





-Q I I {E3) dvdx = -q I / /^s^s dvdx = q j / fv3d-,^(j) dvdx 
/ dv dx 



9tp + g div^r / f „ dv > dx 



B 



j dt{eodiY,E) dx-q [ [ V.0 • ^-^/ dvdx 
-/ V,0-9i(£o^)dx + g / (E)-^^/da; 



2 dt 



implying 



/ ii^;Pdx 

dt I 2 




By Remark 15.21 we know that 

X (Qfpl) if, f) dvdx 






m / dvdx f = 0- 

3 2 



V 

x-\ ) (Qfpl) (/, /) dvdx 



— {Qfpl) (/, /) dt;dx = 
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and therefore, multiplying f l52|) by x one gets 



d f f _ f f (^E) 

' ' fx dvdx — / / / dvdx = 0. 



dt Jm.3Jm.3 JmsJirs -B 

Appealing to the Poisson equation we have 



q f{E) dvdx = q fE dvdx = I pE = Eq I div^EE dx 



= €o j ^ (^div,(E ®E)- ]^V^\E\^^ dx = 

and therefore ^ /ks/ks dvdx = 0. In particular the mean Larmor circle center is 
left invariant 

d_ j lu^Iu^f dt;dx | 
dt I J^J^.f dvdx J 

By Remark 15.21 we know that 

/ / {lo^^\x\^ + 2u:,x ■ ^v) (Qfpl) (/, /) dvdx = 
and for any iIj{x) G we can write 

tp{x) diYx f fvdvdx = —[ [ f V ■ Vx'ip dvdx 

f Tip dvdx = j if) Tf dvdx = 

• Jr3J^3 

saying that div^ f^^ fv dv = 0. Therefore we deduce 

^ /" / 1"''! I 3; . — ^\ j dvdx = [ [ f ^ '^ ■ (x H ^ df dx 




dt Jr3 Jk3\ 2 ojc J jRiJ^s B \ ^ 

I I f—fz- ■ (x H J df dx = — [ I f— ■ """x df dx + f ^ div^ [ fv dv dx 

Jr?J^3 B \ uJc J JrsJks B J^i Buc Jk3 

/ ^-^Ediv^Edx = / ^x- (div^(E® ^) - ^V^l^n dx 
qB J^3 qB Jk3 \ 2 / 

£^ 4. (EiE,) - xi ^ d,^ {E,Ej) I dx 





implying that the mean Larmor circle power (with respect to the origin) is left invariant. 
□ 
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A Proof of Theorem 



1.2 



Proof. The Fokker-Planck kernel being a second order differential operator, we appeal 
twice to Proposition 13.31 For any / G kerT, taking = MV„ (//M) yields 



^V. ( ) ) ) (54) 



M 



div^^3f 





V 


M/ 


\v\ 


^) 


V 


m) 


\v\ 



Since dy^ commutes with (■) (cf. Proposition 13. 2p we deduce 



It remains to compute the averages 



These averages come easily, thanks to Proposition 13. 3[ observing that 



^)-d.v.a,0,0), ^J^)-d.v.(o,^, 



Virl XI ■ = div^ f , V^fXl ■U=div^fXu^ -2 ^ 



We obtain 



Our conclusion follows by combining ( 15^ . ( ISSj) . (!56|) . (l57|l . The diffusion matrix £ is 
positive and for any ^ = (^j,, ^t,) G we have 

with equality iff = = {0, 0) and ^^,3 = 0. 1^ 
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B Proofs of Propositions 15.31 



5.6 



Proof, (of Proposition 15. 3p We follow the same arguments as those in the proof of 



Proposition 14.21 The details are left to the reader. 
1. Using cylindrical coordinates we obtain 

f^2'K P /•TT 



J := / / D{v,v')f{x,v') dv'\ = ^ [1 r [ Di\v\e'^,v,,r'e'^^+'^\v',) 

\Jr3 I ^'^ Jo JR J ~iT JM.+ 

- - i^^^,X3,rV(^+"),t;^) r'dr'd^dt;^da 



X f [x + 

1 p27r p PIT p 

Jo Jm. J -TV Jr+ 

( -'-||u|e*"j -'-/^'e«(¥'+") j \ 
X 5( X H — + — xs, r', ^3 r'dr'dv^dwgda (58) 

where in the last equality we have used the constraint / G kerT i.e., there is g such 
that 

f{x,v) = g (w+-^,X3, \v\,V3^ , {x,v) G X 

We have r'e*^'^"''"-* — |tJ|e*" = le^^'^'^"^ where = + (r')^ — 2rr'cos(.p, r = \v\ and 
(r')2 = + /2 _^ 2rl cos tfj. Notice that G (0,7r) if v? G (0, vr) and ip G (-7r,0) if 
(f G (— 7r,0). Also ip = ipi^ip) is odd with respect to ip that is i'i—ip) = —ip^f). By 
hypothesis we deduce that 

D{re'°',V3,r'e''^'^+"\v's) = D{r,V3,r' ,v'^,<f) re''' + D'{r,V3,r',v's,<f) rV(^+") 

and thus, if we denote by R{a) the rotation of angle a in one gets 

/ / / Ricy + nDir,V3,r',v',,ip)re-'^ + D'ir,V3,r\v'„ip)r'e'^-^-'^^] 

Jo Jk J -it Jr+ 

X (7 ( X H 1 , X3, r', t>o r'dr'do^dtioda 



1 



'^27r /• /'TT 

i?(a) [Z)(r, ^3, r', w^, y;) re"^'^ + Z)'(r, ^^3, ^3, V') ^'e^^"^" 




'0 J-TT 



+ 



X g I X -\ 1 , X3, r', t>3 ) r'dr'dyjdfgda. 
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Using the symmetry of ip with respect to (f and changing 99 against I yield 
Jo Jo Jr+ 

+ D\r, V3, r', ^) rV^'^"'^) + b'{r, V3, r', v'^, -ip) r'e-'^^-"^^] 
X I X H 1 , X3, r', t>3 J r'dr'dv^dfgdo; 

/ / / / R{a)[b{r,V3,r',v'„y:>)re-'^ + b{r,V3,r',v'„-^)re'^ 

'0 7|r-r'| 7lR+ 



TT 



+ D'(r, t;3, r', v'^, ip) rV^'^"^) + D'(r, V3, r', v'^, -ip) r'e''^^^-^^] 
f ^|/e^"| , A r'dr'ldldv'da 



(r-r')2 v/(r + r')2-Z2 
Notice that = ei ® e~*" — 62 ® """e"*" and for any a' G [0, 27r) we have 

x + — + ^ i,X3,r',t;U =/ U+ — + ^ ^-^ i,X3,rV",t;^ 

Performing the change of coordinates v' = (r'e^^'jfg) and —z = ^{/e*"} leads to 



'0 



27r2 

+ D'{<p) rV(^-^) + /)'(-<^) r'e-*(^- 



X / x + — + ^ i - ^ i,X3,rV" x{r,r', -^{le'"})r'dr'da'dv'ldlda 

V UJc J 

= / P(|t;|,t;3, |w'|,t;3,z)/ (x + — - — -,X3,v'] dv'dz 

Jr2 Jk3 V UJc J 



The statement in 2. follows similarly. 1^ 
Proof, (of Proposition I5.6p Observe that 

{f).s = ([ f{x,v')a{\v-v'\)dv') I-lj f{xW)a{\v-v'\) ^''-'''^®%"''^ dv) 

Applying Proposition 14. 21 with C = 1 one gets (Jj^g jo dv'^ = /]k2/r3 di^'dx^dxg. 
It remains to compute the second average. Using cylindrical coordinates and the con- 
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straint Tf = i.e., f{x,v) = g (^x + -^^,2:3, |f 1,^3^ we obtain 

\^M3 \V — V \ 



1 /"^Y r f a( I (|t;|e*"-rV(^+"),t;3-i;^; 



27r Jo Jr+ I ( I e^" - r'e'(^+^) ,V3-v's)\^ 



'V 



xf{x + ^ — i,a;3,rV(^+"\t;^ ) r'dr'dipdv'^da 



X X H — — - + — X3, r', r'dr'dipdv'Aa. 

\ <^c <^c J 

We introduce ^,-0 such that r'e^^'^'^°'^ — \v\e^°' — le^^'^^°'\ We have the relations = 
r^ + (7"')^ — 2rr'cos</7, r = \v\ and (r')^ = r^ + Z^ + 2rZcos'0- Since l,ip are not depending 
on a one gets 



X gf I a; H 1 ,a;3,r',t'o ) r'dr'df^dvodo! 



27r io Jm i-^ iR+ + (^^3 - v'^y 
y. g[x A 1 , X3, r', J r'dr'dt/^di'gdQ;. 

\ UJc OJc J 

Since Z is even with respect to </? we obtain, after changing (/? against I 
X g \ x -\ 1 , X3, r', vL ) r'dr'dc/^dvqdo; 

Jin f-^l .(VTO^) 

TT io J^Jlr-r'l Jr+ ^ + (^3 " ^'s) 

/ ^{/e'"} , A r'dr' Idl dv'da 

For any a' e [0, 27r) we have 

y x+ — + ^ ^-,X3,r',v'A = f {x+ — + -^ ^---^ i,X3,rV",^;^ 
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Performing the change of coordinates v' = (r'e*" , Wg) and z = — ^{/e*"} leads to 

"2-71/' r r^TT , 



K 



'0 Jr+JrJo Jr+ ' +11^3 t^sJ 




X / ( x + — + ^ - ^^^^ — ^,X3,rV" x(^,^',-^{^e'"})r'dr'da'dt;^MMa 

1 12 , / TTi ^' ^3 - ^3 / U + , X3,v ]x dv dz 

+(^3-^3) V (^c tJc J 

(VkP + (t;3-t;02)/(x;,x'2,X3,t;') X r^^T^T ^ dt;'da;;d4. 





Lol \ la 



Finally we obtain 




{f)as = ^c / o^(VkP + (^^3 -^^3)^ )fix[,X2,X3,v') X Si (^z, t;3 - f 3) ) dv'dx'^dx'^. 



□ 



C Proofs of Propositions 15. 7L 15.8 



Proof, (of Proposition I5.7p Let us introduce the notation 

Ux,v)= f a{\v-v'\)S{v-v')f{x,v')VJ-{x,v)dv'. 
Thanks to Proposition 13.31 we have 

(Q^rPLifJ)) = (div. Q = ^div,{(^e.> + ^) ^ - ^) ^} 

We need to compute {^y) , {C,v ■ v) , {C,v ■ ^v). By Proposition l3. ll we know that Yl^=o 
^ x,vipi = I and thus 

5 



Similarly we have 

5 



i=0 
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leading to 



Uc\v\^ \ Uc J \v\ 

Taking into account that all derivations 6* ■ Vxv,^ < ^ < 5 leave invariant kerT, cf. 



Proposition 13. we obtain 

(e.) = (/, ( ^v, 0))^, + {fXs (-^-^m) + if, {V, 0)), ■ ^-^^ 



(/'(^^'0)>.5^ A^,0) (^T;,0)\ (/,(tJ,0))^,^ ^(^^,0) (^,0) 



1^1 \ |f| |f| / \v\ \ \v\ \v\ 

where the lines of the matrix E G A^3(M) are 62,— ei,0. Similarly, thanks to the 
identities (/, (v, 0), ( -^v, 0))^^ = (/, ( -^v, 0), (tJ, 0))^^ = we obtain 



v\ / \ \v\ 



V\ \ COr / \V 



if, (v, 0), {v, 0))^, ( ( _ ^ /^ (/, ^ CMMk 



and 



\v\ 



I^P Wclfl \v\ \ Uc ' 

(/,(^U,0),(^U,0)>^ /(l;,0) (^l;,0)\ , /^^(/,(^^,0)>^^ 



v\^ \ \v\ ' 



+ E^^^^—^^, 0]}- V.,.,./. 



In the last equality we have taken into account that 

es ® ^ ,_ ^^"^ = 0. 
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Clearly {Qpp^if, /)) has the form in (139|) with = | where 



A+ A 



vv 



\ xxy xv) 



\Z\ \Z\ \ \v\ ' \v\ 



{r — r' cosip){v3 — v'^y {z,0) / ( , 0) {v,0) 



+ 1 



2;| [ I^P + (^3 — fg)^ ] \z\ \ \v\ ' \v\ 



+ (f3 — fg)^/ |f| \ |f| \v\ 

r — r' cos (/9 {v, 0) ^ f {^z, 0) 



(r -r'cosy?)^ ^^ (^I;,0) /^(^t^,0) (v,0) 



+ ("^3 ~ "^s)^/ l"^! V l"^! 

+ i^z, v', - vs) )E, ES{ i^z, v', - vs) )e3 ® 63) 



and 



r — r' cos 09 (^v,0) f (^z,0) 
pI |f| \ \z\ 

_ ^ jr-r'cosipy \ (U,0) ^ ^(^I;,0) (U, 0) 



^ ( r — r cos 



^){v3-v',) {v,0) ((^3-t^3)l|^,-k|e3) 



zW\z\^ + {v3-v',y \V\ ^|^|2 + (^3_^/)^ 

_^ (r - /cosv?)(t;3 - t;^) ^ Z' ( "^^^ 0) _ (M) 

|zp + (f3 — v'^y \ \v\ ' \v\ 

+ (-63 ® e35( {^z, - V3) )E, 63 ® e35( (^z, - V3) )e3 ® 63). 
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It is easily seen that the matrix A'^ writes 



+ 



(r')^ sin^ (/? 
+ {vs - v'^y 
(r — r' cos (/?)^ 

r — r' cosip f {^z,0) 
~\ i ' 



iv,0) (%,0) 



-v,0) (y,0) 



02 



v\ \v\ 



+ 



r — r' cos / (f , 0) (^f,0)\ /( 



\v\ 

\z\ 



,0 



(r -r'cosv9)(v3 -Vs) ('^^s -t'3)l^>-k|e3) /(^i;,0) (i;, 0) 



(r - r' cos if) (vs-v's) ({ ^v, 0) {v, 0) \ ((^3 - ^^3) 



\V 



(61) 



where 



B- 



'ES{ i^z, v', - vs) )E ES{ i^z, v', - V,) )es ® 63 

-63 (8) e35'( {-^z, - V3) )£; 63 (g) e35'( (-^z, - V3) )e3 (8) 63 

Observe that for any 2; G M^, 1)3, G M the family 



iz,0) 



-z, v'^ - V3) 



y|^|2 + (^3 _ ^/)2 ' ^|^|2 + (^3 _ ^/)2y ' ^|^|2 + (^3_^/)2 

is a orthonormal basis of M^. Therefore we have 



S{ ( z. Ug - i;3) ) 



- f g 



^|^|2+(^3_^/)2|^|'^|^|2 + (^3_^/)2 



implying that for any i = {ix,iv) 



-z,0) 



,0 + 



^3 -^-s 



[z,0) 



k|e3 



02 



^|;,|2 + (^3_^/)2 ' ^|^|2 + (^3_^/)^ 



and thus 



B+^ 



-z,0) 



,0 



fg -fg 



{Z,0) 



2 U; 









2 + (y^-y'^)2 



(62) 
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Observe that 



At + A+ + At + B+ 



{r'y sin' if {V3 - v'^y f{v,0) ( ^U, 0) 



+ 



zP[ l^;^ + (^3 — t;3)2 ] \ \v 

r — r' cos ip f (v, 0) ( ^v, 0) 

\v\ ' \v\ 



+ 



-z,0) 



since (l^lifip)^ = 1 _ ilDl^ and 



At + At + A+ + B' 



(rO^sinV A ^,0) (tJ, 0) 



(r-r'cos(p)(t;3-^^3) A^^.O) (^.0)\ ((^'s - I'a)!^. -k^a 



+ 



Our conclusion follows by combining (16T|) . ( 162|) . (163|) . ( 16^ . 



Proof, (of Proposition I5.8p We consider 

^^{x,v)= / a{\v-v'\)S{v-v')f{x,v)V^>f{x,v')dv'. 



By Proposition 13.31 we have 



(Q^P.(/,/)> = (div. Q = ^div,{(^e.> + {^--^)^\- (^-- ^) ^} 
f / ^tJX '^v I u \ u 1 



J. [/. ( ^t;, 0) \ ( ^tJ, 0) I (U,0)\(tJ,0) , 

As in the proof of Proposition 15.71 we obtain 

V„'/ X, t; = 6^ ■ V^y/ — + ,d^' f ] + b ■ V^y/ —=- 

Uc\v'\' \ Uc ^ J \v'\ 
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and therefore 



+ fix, v) (9,, /, e3>^^ + fix, v) C—^^, 0) 



<tS 

2 I I r, X r' - r COS 0) / (w', 0) (^'y',0)\ ^j/j/j/ 



^-^c / n rr ® ^t^' — — ^ujcxyf dv'dx[dx2 

Jr^Jr^ \z\ \z\ V 1^ I 1^1 

^ — , ^ojcx,v'f dv dx^dx2 

\v'\ \v'\ ) 

-'A \ \ <^f{x,v)xS{ (^^,^^3 - ^3) ){E, -63 ® e3)Va;,xy/ d^;'da;^da;2• 
Similarly, thanks to the identities 



/|2 



<0),(l;,0)\ =/^-^^,KO),(^t;,0)\ 

cl-l /aS \ Pl /aS 



we obtain 



■ ) = V) { — — , {v, 0) ) -61 + /(x, v) { — — , [v, 0) ) -62 



+ /(a;, -y) ( , [v, 0) ) • 63 + /(x, -y) ( —= , {v', 0), {v, 0) 



/ (tS \ 1^11^ I / aS 




= -^c / _ / '^/^ + (^3-^^)2 ^' J • V.,.y/ d^ dx,dX2 

9 / / „ / (r — r'cos</?)(r' — rcos(i3)\ / ( -"-v', 0) (v', 0) 



^iOox,v'f dw'da;ida;2 



and 



J,, ./d^J{x,v') .x_„A , J,, ./d,^J{x,v') ^_ 



+ /(a;, -y) ( , ( 0) ) ■ 63 + /(x, i;) ( - ,( ^;^0),( i;, 0) 



^1 / aS \ <^c\v\\V\ / 

-^l I I 'yfx '""Y^~" ) • V...,,/ dv'dx\dx', 

+ / / <^/X [coscp- ^ ,^ • V^.^y/ d^;'dx;dx^. 




(^^3 --^3)^ 
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A~ A~ 

Obviously {Qppiifj f)) has the form in (HOj) with A~ = j xv \ ^y^qiq 



A- A- 

vx vv 



r;-rcos<^ (^z,0) _ /KO) (V,0) 



{r' -rco^^){v^-v'^y {z, 0) / ( ^v', 0) (t;', 0) 
~r ~, — rn — it; 7 rm~ — i — i — 'o) 



z\[\z\'^ + {v^ — v'^Y ] \z\ \ \v'\ ' \v'\ 

rr' sin^ \ (IJ, 0) / 0) ( ^t7, 0) 



+ COS(y9 



+ ("^3 — "^s)^/ 1"^ 



r — r' COS 09 (f , 0) / (-'-2;, 0) „ 

+ n — - ^T=r- ® I I , 

\z\ \v\ \ \z\ 

( (r — r' COS (r' — r COS \ ( -'-U, 0) / ( 0) (f',0)' 

+ + — I 12 I ^ Mi — i-i ® ^ — , — 

_ {r -r' co^if>){v'i-v'^ ( -^U, 0) ^ / ^3 - ^3 (^^ 0) _^ 

+ ("^^3 ~ "^^3)^ 1^1 V \A \A 

- ES{ {^z, v', - V,) ){E, -eg ® eg) (65) 



and 



_ / ^tW^\ ( ^U, 0) ^ f {v', 0) ( ^t;', 0) 

[A^,, A^J - ^cos ^ - ^^—^-—^ J ^ 

r — r' cos 09 (^v,0) / (^z,0) „ 

+ r-< — - I- ® II , 

|z| l^;! \ pI 

(r — r' cos V9) (r' — r cos </)) \ (t^, 0) f(^v',0) (.v',0) 
cos V9 H — 1^1 



+ (^3 - ^3)^ 

(r -r'cosv9)(t^3 -^^3) (^^ 0) ^ /^^3 - ^^3 (^^ 0) 
1 12 1 ^ M2 PF — n TF' 

1^1 + 1^3 ~ ^3) V 1-2^1 1-2^1 

(r' — r COS (y9)(t>3 — t>3) / ( -"-f', 0) (f 0) 

63 ^ 



+ (^^3 - ^^3 



^2 



-eg® eg S{ {^z,v'3-V3) )(^,-eg®eg). (66) 
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It is easily seen that the matrix A writes 

U,0) (^^7,0) 



A 



C0SV9 



rr' sin^ ^ 



{v',0) (V,0) 



\v'\ 



(r — r' cos (p)(r' — r cos ip)\ f (^v,0) (v,0) 

r' -rcos(^ /(^z,0) \ ( V, 0) 

~; — ; — (S) — F=; — , — — 



-v',0) (f',0) 



+ 



\z\ \ \z 
r — r' cos if ( (l;, 0) ( ^v^ 0) 



\v\ ' \v\ 



\V'\ \v 



-2 



+ 



(r^ - rcos(^)(t;3 - v'^) (,(^3 ~ ^3)l^> 'W?,) / ( V, Q) _(<0) 



5- 



lea 



^|^|2+(^3_^/)2 



\v'\ 



(67) 



where, cf. 



-63 ® 635* ( (-^2;, - ^73) )E 63 (g) e35'( (-^z, 1^3 - ^73) )e3 ® 63 



-z,0) 







^3 -^^3 



:^,o) 



\z\e3 



^\Z\^ + {V,-V',y \Z\ Vkl' + (^3-^3)' 



Observe that 



+ 



rr' sin^ (/9 (f3 — f3)^ f{v,0) {^v,0) 

|z|2 + (t)3 — t>^)2 j \^ 1^1 ' 1^1 

r — r' cos(p f (v,0) {^v,0)\ f {^z,0) 



[v',0) (V,0) 



\v'\ 



\v\ ' Ifl 



+ 



r cos if — r' f (t>', 0) ( -"-f ', 0) 

I — 7 1 : I — Ti I 



z,0) 



(69) 



since 



rr' sin <^ (^3 — f 3) ^ (r — r' cos y?) (r cos — r') 



rr' sin^ 



|z|2 + (t)3 — t)02 j 



A"^ + (^3 - ^3)' 
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and 



A:; + A7 + A 



6 ^ -^2 



rr' sin V /( 0) {v,Q)\ [{^v\Q) (i;',0) 



+ 



(r — r' cos V9)(i;3 — fg) / ( -'-f , 0) ("y, 0) 



If I 



(^3-^)^,-k|e3 



1^3 



(rcos(/?-r')(f3-f3) /( V,0) (f',0)\ ((^3 - ^^3)^^ 

klVkP + (^3^^ V PI ' W) ^\z? + {v^-vW 



(70) 



since 



rr' s\.v? if {r' cos if — r){r' — r cos if) iy 2, — v'.^Y {r — r' cosip){r' — r cosip) 

[ + ("^3 — "i^s)^ ] ^ + ("^3 — I's)^ 

Our conclusion follows by combining (167|) . ( 168|) . (!69|) . ( ITOj) . ^ 
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